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1 Introduction 



Strong interactions as described by Quantum Chromo dynamics (QCD) is a remarkable branch 
of physics where the observable entities - hadrons and nuclei - are very far from quarks and 
gluons in terms of which the theory is formulated. To make matters worse, the scale of strong 
interactions 1 fm is nowhere to be found in the QCD Lagrangian. If we restrict ourselves to 
hadrons 'made of u, d, s quarks and glue, the masses of those quarks can be to a good accuracy 
set to zero. In this so-called chiral limit the nucleon is just 5% lighter than in reality. In the 
chiral limit there is not a single dimensional parameter in the QCD Lagrangian. The 1 fm 
scale surfaces via a mechanism named the 'transmutation of dimensions' . QCD is a quantum 
field theory and beeing such it is not defined without introducing of some kind of ultra-violet 
cutoff /i. There is also a dimensionless gauge coupling constant given at that cutoff Q;s(/i)- The 
dimensionful quantity A determining the scale of the strong interactions is the combination of 
// and Q;s(/i): 



h = -N,--Nf, b2 = -N'--N,Nf + ^, (1.2) 
1 3 3-^' 3 3 A^c 

where Ac = 3 is the number of quark colours and Nf is the number of acting quark fiavours. 
The ultra-violet cutoff n sets in the dimension of mass but the exponentially small factor makes 
A much less than fi. To ensure that A is actually independent of the cutoff, one has to add 
that as{iJ>) has to decrease with /i according to 

2- ,.,„ii + J^„.„^ + om. (1.3) 



asin) A 2fei A2 \\nf^ 

This formula is called 'asymptotic freedom': at large scales ag decreases. 

All physical observables in strong interactions, like the nucleon mass, the pion decay constant 
F^, total cross sections, etc. are proportioal to A in the appropriate power. That is how the 
strong interactions scale, 1 fm, appears in the theory. One of the theory's goals is to get, say, 
the nucleon mass in the form of eq. (1.1) and to find the numerical proportionality coefficient. 
Doing lattice simulations the first thing one needs to check is whether an observable scales with 
as as prescribed by eq. (1.1). If it does not, the continuum limit is not achieved. In analytical 
approaches, getting an observable in the form of eq. (1.1) is extremely difficult. It implies doing 
non-perturbative physics. The only analytical approach to QCD I know of where one indeed 
gets observables through the transmutation of dimensions is the approach based on instantons, 
and it will be the subject of this paper. 

Instantons are certain large non-perturbative fluctuations of the gluon field discovered by 
Belavin, Polyakov, Schwartz and Tyupkin in 1975 [1, 2], and the name has been suggested in 
1976 by 't Hooft [3], who made a major contribution to the investigation of the instantons 
properties. The QCD instanton vacuum has been studied starting from the pioneering works 
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in the end of the seventies [4, 5]; a quantitative treatment of the instanton ensemble has been 
developed in refs. [6, 7]. The basic ideas of the instanton vacuum are presented in section 2. 

Instantons are not the only possible large non-perturbative fluctuations of the gluon fleld: 
one can think also of merons, monopoles, vortices, etc. I briefly review that in section 3 
where also certain new material on dyons is presented. However, instantons are the best 
studied non-perturbative effects. It may happen that they are not the whole truth but they are 
dcflnitcly present in the QCD vacuum, and they are working quite effectively in reproducing 
many remarkable features of the strong interactions. For example, instantons lead to the 
formation of the gluon condensate [8] and of the so-called topological susceptibility needed to 
cure the U (1) paradox [3, 9]. The most striking success of instantons is their capacity to provide 
a beautiful microscopic mechanism of the spontaneous chiral symmetry breaking [10, 11, 12, 13]. 
Moreover, instantons enable one to understand it from different angles and using different 
mathematical formalisms. These topics are central in the review and are presented in sections 
4,5 and 6. 

We know that, were the chiral symmetry of QCD unbroken, the hghtest hadrons would 
appear in parity doublets. The large actual splitting between, say, , 1535) and N{^'^, 940) 
implies that chiral symmetry is spontaneously broken as characterized by the nonzero quark 
condensate < ipip >~ — (250MeV)^. Equivalently, it means that nearly massless ('current') 
quarks obtain a sizable non-slash term in the propagator, called the dynamical or constituent 
mass M(p), with M(0) ~ 350 MeV. The p-meson has roughly twice and nucleon thrice this 
mass, i.e. are relatively loosely bound. The pion is a (pseudo) Goldstone boson and is very light. 
The hadron size is typically ~ 1/M(0) whereas the size of constituent quarks is given by the 
slope of M{p). In the instanton approach the former is much larger than the latter. It explains, 
at least on the qualitative level, why constituent quark models are so phenomenologically 
successful. 

It should be stressed that literally speaking instantons do not lead to confinement, although 
they induce a growing potential for heavy quarks at intermediate separations [14]; asymptoti- 
cally it flattens out [4]. However, it has been realized long ago [5, 15], that it is chiral symmetry 
breaking and not confinement that determines the basic characteristics of nucleons and pions 
as well as their first excitations. After all, 99% of the mass around us is due to the spontaneous 
generation of the quark constituent mass. Probably one would need an explicit confinement 
to get the properties of short-living highly excited hadrons. According to a popular wisdom, 
moving a quark away from a diquark system in a baryon generates a string, also called a flux 
tube, whose energy rises linearly with the separation. This is expected in the "pure-glue" world 
with no dynamical quarks. However, in the real world with light quarks and the spontaneous 
chiral symmetry breaking the string energy exceeds the pion mass = 140 MeV at a modest 
separation of about 0.26 fm, see Fig. 1. At larger separations the would-be linear potential is 
screened since it is energetically favourable to tear the string and produce a pion. Virtually, 
the linear potential can stretch to as much as 0.4 fm where its energy exceeds 2171-^ but that 
can happen only for a short time of 1/m^. Meanwhile, the ground-state baryons are stable, 
and their sizes are about 1 fm. The pion-nucleon coupling is huge, and there seems to be no 
suppression of the string breaking by pions. The paradox is that the linear potential of the 
pure glue world, important as it might be to explain why quarks are not observed as a matter 
of principle, can hardly have a direct impact on the properties of lightest hadrons. 
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Figure 1: The lattice-simulated potential between static quarks in pure glue theory [16] exceeds 
at the separation of 0.26 fm (left). The screening of the linear potential by dynamical quarks 
is clearly seen in simulations at high temperatures but below the phase transition [17] (right). 
As one lowers the pion mass the string breaking happens at smaller distances; the scale is 
v/a~ 425MeV~ (0.47fm)-i. 

Even for highly excited hadrons lying on linear Regge trajectories the situation is not al- 
together clear. The usual explanation of resonances lying on linear trajectories is that they 
are rotating confining flux tubes attached to quarks at end points moving with the speed of 
light. Why should 350 MeV quarks bound by a ^Ja ~ 425 MeV string move with the speed 
of light is not clear, but if they do not, the trajectories are not linear. In this picture, the 
finite (and experimentally large) width of resonances is due to the same string breaking by 
meson production. However, if there is no string in the ground-state nucleon, why should it 
be excited in a collision? The lightest degrees of freedom in the real world are pions and one 
might expect that they are the first to be excited. An alternative explanation of resonances 
lying on linear trajectories is that they are rotating elongated lumps of pion field [18], and their 
decay is due to the normal pion radiation. It follows then that the dominant decay of a large-J 
baryon resonance is a cascade Barj Barj_i + tt ^ Barj_2 + vrvr — > ... whereas if it is due 
to the string breaking it rather has a different pattern Barj Bar^j/2 + Mes^j/2. Studying 
the decay patterns of high-J resonances could be illuminating for understanding the relation 
between confining and chiral forces. 

Leaving aside the unsettled question of highly excited resonances, the situation with the 
lightest and most important hadrons 7r,p, A^, A... is, to my mind, clear: it is the spontaneous 
chiral symmetry breaking (SCSB) rather than the expected linear confining potential of the pure 
glue world which is the key to understanding of their properties. Therefore, since the instanton 
vacuum describes successfully the physics of the chiral symmetry breaking, one would expect 
that instantons do explain the properties of light hadrons, both mesons and baryons. Indeed, a 
detailed numerical study of dozens of correlation functions with different quantum numbers in 
the instanton medium undertaken by Shuryak, Verbaarschot and Schafer [19] (earlier certain 
correlation functions were computed analytically in refs. [11, 12]) demonstrated an impressing 
agreement with the phenomenology [20] and with direct lattice measurements [21], see ref. [22] 
for a review. In fact, instantons induce strong interactions between quarks, leading to bound- 
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state baryons with calculable and reasonable properties. There are specialized reviews on this 
subject, therefore I touch it only briefly here (sections 8 and 9). 

More recently, there has been much activity in applying instantons to explain various phe- 
nomena in high energy processes including heavy ion collisions. For that reason, 1 have included 
section 7 which suggests a new point of view on the pomeron which might be also related to 
instantons. 



2 What are instantons? 

2.1 Periodicity of the Yang-Mills potential energy 

Being a quantum field theory, QCD deals with the fiuctuating gluon and quark fields. A 
fundamental fact [23, 24] is that the potential energy of the gluon field is a periodic function 
in one particular direction in the infinite-dimensional functional space; in all other directions 
the potential energy is oscillator-hke. This is illustrated in Fig. 2. 

To observe this periodicity, let us temporarily work in the — gauge, called Weyl or 
Hamiltonian gauge, and forget about fermions for a while. The remaining pure Yang-Mills or 
"pure glue" theory is nonetheless non-trivial, since gluons are self-interacting. For simplicity I 
start from the SU (2) gauge group. 

The spatial YM potentials (x, t) can be considered as an infinite set of the coordinates 
of the system, where i = 1,2,3, a = 1,2,3 and x are "labels" denoting various coordinates. 
The YM action is 

s = ^ jd'^ = / * - 1^ B'^) (2. 1) 

where E is the electric field stregth, 

£;»(x,i)=i»(x,i) (2.2) 
(the dot stands for the time derivative), and B is the magnetic field stregth, 

Bti^, t) = \e^,k {d,Al - d^A] + e^'^A)Ai) . (2.3) 

Apparently, the first term in cq. (2.1) is the kinetic energy of the system of coordinates 
{A^(x, t)} while the second term is minus the potential energy being just the magnetic energy of 
the field. The simple and transparent form of cq. (2.2) is the advantage of the Weyl gauge. Upon 
quantization the electric field is replaced by the variational derivative, -Ef (x) —ig^b jbA^ix)^ 
if one uses the 'coordinate representation' for the wave functional. The functional Schrodinger 
equation for the wave functional ^[^"(x)] takes the form 
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where 8 is the eigenenergy of the state in question. The YM vacuum is the ground state of the 
Hamiltonian (2.4), corresponding to the lowest energy £. 

Let us introduce an important quantity called the Pontryagin index or the four-dimensional 
topological charge of the YM fields: 

1 r ~ ~ 1 

Qt = j d'^x F^^F^^, F^^ = -e^uapF^fj- (2.5) 

The integrand in eq. (2.5) happens to be a full derivative of the four-vector X^: 

^F;,F;, = d,K,, K, = [AldpA^^ + ^e«^^^«4A^) . (2.6) 

Therefore, assuming the fields are decreasing rapidly enough at spatial infinity, one can 
rewrite the 4-dimensional topological charge (2.5) as 

Qt = / d'x(9oXo - diK,) = j dtjj d^xKo. (2.7) 
Introducing the Chern-Simons number 

Ncs = / (^'x Ko = J d^x e'^^ (Ald^Al + ^e^^^A^AjA^) (2.8) 

we see from eq. (2.7) that Qt can be rewritten as the difference of the Chern-Simons numbers 
characterizing the fields &X, t — ±oo: 

Qt^Ncs{+<^)-Ncs{-<^). (2.9) 

The Chern-Simons number of the field has an important property that it can change by 
integers under large gauge transformations. Indeed, under a general time-independent gauge 
transformation , 

Ai^U^AiU + iU^diU, A = A1^, (2.10) 
the Chern-Simons number transforms as follows: 

Ncs ^Ncs + Nw + ^j d^x e'^^d, Tr {d^VV'^A^). (2.11) 

The last term is a full derivative and can be omitted if, e.g.^ Ai decreases sufficiently fast at 
spatial infinity. is the winding number of the gauge transformation (2.10): 

Nw^^J d'x e'i^ [{U^d,U){U^d,U){U^dkU)\ . (2.12) 

The SU{2) unitary matrix U of the gauge transformation (2.10) can be viewed as a map- 
ping from the 3-dimensional space onto the 3-dimensional sphere of parameters 5"^. If at spatial 
infinity we wish to have the same matrix U independently of the way we approach the infin- 
ity (and this is what is usually assumed), then the spatial infinity is in fact one point, so 
the mapping is topologically equivalent to that from to S^. This mapping is known to be 
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Figure 2: Potential energy of the gluon field is periodic in one direction and oscillator-like in 
all other directions in functional space. 



non-trivial, meaning that mappings with different winding numbers are irreducible by smooth 
transformations to one another. The winding number of the gauge transformation is, analyt- 
ically, given by eq. (2.12). As it is common for topological characteristics, the integrand in 
(2.12) is in fact a full derivative. For example, if we take the matrix f/(x) in a "hedgehog" 
form, U = exp[i{r ■ T)/r P{r)], eq. (2.12) can be rewritten as 



Nw = - dr— sm^ P = - 



^ sin2P 



. integer (2.13) 

n J dr TT L 2 J 

since P(r) both at zero and at infinity needs to be multiples of vr if we wish U{f) to be 
unambigiously defined at the origin and at the infinity. 

Let us return now to the potential energy of the YM fields, 

V=^Jd^^iBt)\ (2.14) 

One can imagine plotting the potential energy surfaces over the Hilbert space of the co- 
ordinates A'}{x). It will be some complicated mountain country. If the field happens to be 
a pure gauge, Ai = iWdiU, the potential energy at such points of the Hilbert space is natu- 
rally zero. Imagine that we move along the "generalized coordinate" being the Chern-Simons 
number (2.8), fixing all other coordinates whatever they are. Let us take some point A'^{x.) 
with the potential energy V. If we move to another point which is a gauge transformation of 
A^{x.) with a winding number Nw, its potential energy will be exactly the same as it is strictly 
gauge invariant. However the Chern-Simons "coordinate" of the new point will be shifted by 
an integer Nw from the original one. We arrive to the conclusion first pointed out by Faddeev 
[23] and Jackiw and Rebbi [24] in 1976, that the potential energy of the YM fields is periodic 
in the particular coordinate called the Chern-Simons number. 

2.2 Instantons in simple words 

In perturbation theory one deals with zero-point quantum-mechanical fluctuations of the YM 
fields near one of the minima, say, at Ncs = . The non-linearity of the YM theory is taken 
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into account as a perturbation, and results in series in where g is the gauge couphng. In that 
approach one is apparently missing a possibility for the system to tunnel to another minimum, 
say, at Nqs = 1 • The tunneling is a typical non-pcrturbative effect in the coupling constant. 

Instanton is a large fluctuation of the gluon field in imaginary (or Euclidean) time corre- 
sponding to quantum tunneling from one minimum of the potential energy to the neighbour one. 
Mathematically, it was discovered by Belavin, Polyakov, Schwarz and Tyupkin; [1] the tunnel- 
ing interpretation was given by V. Gribov, sec [2]. The name 'instanton' has been introduced 
by 't Hooft [3] who studied many of the key properties of those fiuctuations. Anti-instantons 
are similar fluctuations but tunneling in the opposite direction in Fig. 2. Physically, one can 
think of instantons in two ways: on the one hand it is a tunneling process occuring in time, on 
the other hand it is a localized pseudoparticle in the Euclidean space. 

Following the WKB approximation, the tunneling amplitude can be estimated as exp(— S*), 
where S is the action along the classical trajectory in imaginary time, leading from the minimum 
at Ncs = at t = — oo to that at Ncs = 1 at t = +oo. According to eq. (2.9) the 4-dimensional 
topological charge of such trajectory is Qt = 1. To find the best tunnehng trajectory having 
the largest amphtude one has thus to minimize the YM action (2.1) provided the topological 
charge (2.5) is fixed to be unity. This can be done using the following trick [1]. Consider the 
inequality 

= J d^x (2F^ - 2FF) = 8/5 - 647r^QT, (2.15) 
hence the action is restricted from below: 

S>%Qr^ (2-16) 

r 9 

Therefore, the minimal action for a trajectory with a unity topological charge is equal to 
871^/51^, which is achieved if the trajectory satisfies the self-duality equation: 

f;^ = f;^- (2.17) 

Notice that any solution of eq. (2.17) is simultaneously a solution of the general YM equation 
of motion D^^F^^ — 0: that is because the "second pair" of the Maxwell equations, D'^^'Fj^^ — 0, 
is satisfied identically. 

Thus, the tunneling amplitude can be estimated as 

A ^ exp(-Action) = exp J c^^^-^ij = exp (^~~^^ = (~^) " ^^--^^^ 

It is non-analytic in the gauge coupling constant and hence instantons are missed in all orders 
of the perturbation theory. However, it is not a reason to ignore tunneling. For example, 
tunneling of electrons from one atom to another in a metal is also a non-perturbative effect but 
we would get nowhere in understanding metals had we ignored it. 
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2.3 Instanton configurations 



To solve eq. (2.17) let us recall a few facts about the Lorentz group 5*0(3, 1). Since we are 
talking about the tunneling fields which can only develop in imaginary time, it means that we 
have to consider the fields in Euclidean space-time, so that the Lorentz group is just 50(4) 
isomorphic to SU(2) x SU{2). The gauge potentials belong to the (|, |) representation of 
the SU{2) X SU{2) group, while the field strength F^i, belongs to the reducible (1,0) + (0, 1) 
representation. In other words it means that one linear combination of F^,^ transforms as a 
vector of the left SU{2), and another combination transforms as a vector of the right SU{2). 
These combinations are 

= VU^,. + K^), Fi = ni{F,. - F,,), (2.19) 

where r], fj arc the so-called 't Hooft symbols described in ref . [3] , see also below. We see 
therefore that a self-dual field strength is a vector of the left SU (2) while its right part is zero. 
Keeping that experience in mind we look for the solution of the self-dual equation in the form 

Ai-n;.x.^-^. (2.20) 

Using the formulae for the rj symbols from ref. [3] one can easily check that the YM action can 
be rewritten as 




(2.21) 



This can be recognized as the action of the double- well potential whose minima lie at $ = ±1, 
and T plays the role of "time" ; p is an arbitrary scale. The trajectory which tunnels from 1 at 
T = — oo to — 1 at r = -|-oo is 

$ = -tanh(^0, (2.22) 

and its action (2.21) is 5 = STr^/gf^, as needed. Substituting the solution (2.22) into (2.20) we 

get 

The correspondent field strength is 



and satisfies the self-duality condition (2.17). 

The anti-instanton corresponding to tunneling in the opposite direction, from Nqs = 1 to 
^cs = , satisfies the antz-self-dual equation, with F —F; its concrete form is given by 
eqs.(2.23, 2.24) with the replacement fj ^ r]. 
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Eqs.(2.23, 2.24) describe the field of the instanton in the singular Lorenz gauge; the singu- 
larity of Afj, at = is a gauge artifact: the gauge- invariant field strength squared is smooth 
at the origin. Formulae for instantons are more compact in the Lorenz gauge, and I shall use 
it further on ^. 

2.4 Instanton collective coordinates 

The instanton field, eq. (2.23), depends on an arbitrary scale parameter p which we shall call the 
instanton size, while the action, being scale invariant, is independent of p. One can obviously 
shift the position of the instanton to an arbitrary 4-point - the action will not change either. 
Finally, one can rotate the instanton field in colour space by constant unitary matrices U. For 
the SU{2) gauge group this rotation is characterized by 3 parameters, e.g. by Euler angles. For 
a general SU{Nc) group the number of parameters is A^^ — 1 (the total number of the SU (iVc) 
generators) minus {N^ — 2)^ (the number of generators which do not affect the left upper 2x2 
corner where the standard SU{2) instanton (2.23) is residing), that is ANc — 5. These degrees 
of freedom are called instanton orientation in colour space. All in all there are 

4 {centre) + 1 (size) + {ANc - 5) (orientafAons) = AN^ (2.25) 

so-called collective coordinates desribing the field of the instanton, of which the action is inde- 
pendent. 

It is convenient to indroduce 2x2 matrices 



aj = (±ia^,l), x^^x^a^, (2.26) 



such that 



2zr'^C = - <^>;, 2iT"C = a;a+ - a-a+, (2.27) 

then the instanton field with arbitrary center z^, size p and colour orientation U in the SU{Nc) 
gauge group can be written as 

2 

~{x- z)2[p2 + {x- zY] ' V— ; - 2 



-ip'^U[a-(x- z)+ - (x- z)JW , 1 , 

A, = A^e = ,!!Ur.2 ^r. LJ ^ Tr (n") = -J^\ (2.28) 



or as 



Al^ ^ ^o^ o . 0"' = Tr(;7tt'^t/a'), O^^'O"" = S""" . (2.29) 

{x - zy[p^ + {x - zyy \ n v / 

This is the explicit expression for the 4A^c-parameter instanton field in the SU {Nc) gauge theory, 
written down in the singular Lorenz gauge. 

^Jackson and Okun [25] recommend to call the d^Afj, = gauge by the name of the Dane Ludvig Lorenz 
and not the Dutchman Hendrik Lorentz who certainly used this gauge too but several decades later. 
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2.5 Gluon condensate 



The QCD perturbation theory imphes that the fields ^"(x) are performing quantum zero-point 
oscillations; in the lowest order these are just plane waves with arbitrary frequences. The 
aggregate energy of these zero-point oscillations, (B^ + E^)/2, is divergent as the fourth power 
of the cutoff frequency, however for any state one has (-F^,^) = 2(B^ — E^) = 0, which is just 
a manifestation of the virial theorem for harmonic oscillators: the average potential energy is 
equal the kinetic one (I am temporarily in the Minkowski space) . One can prove that this is also 
true in any order of the perturbation theory in the coupling constant, provided one does not 
violate the Lorentz symmetry and the renormalization properties of the theory. Meanwhile, we 
know from the QCD sum rules phenomenology that the QCD vacuum posseses what is called 
gluon condensate [8]: 

^ -XFIFI) = -1-(B2 - E^) ~ (200 MeVY > 0. (2.30) 



327r2^ M^'/ 16^2 

Instantons suggest an immediate explanation of this basic property of QCD. Indeed, instanton 
is a tunneling process, it occurs in imaginary time; therefore in Minkowski space one has 
£'f = ±i-Bf (this is actually the duality eq. (2.17)). Therefore, during the tunneling B^ — E^ 
is positive, and one gets a chance to explain the gluon condensate. In Euclidean space the 
electric field is real as well as the magnetic one, and the gluon condensate is just the average 
action density. Let us make a quick estimate of its value. Let the total number of instantons 
and anti-instantons (henceforth 7's and J's for short) in the 4-dimensional volume V be N. 
Assuming that the average separations of instantons are larger than their average sizes (to be 
justified below), we can estimate the total action of the ensemble as the sum of invidual actions 
(see eq. (2.16)): 

{FlJ)V = J d'xFl^ ^ N ■ 327^^ (2.31) 

hence the gluon condensate is directly related to the instanton density in the 4-dimensional 
Euclidean space-time: 

^ -XFIF",^) ^^^ll- (2.32) 



In order to get the phenomenological value of the condensate one needs thus to have the average 
separation between pseudoparticles [8, 5] 

R ~ = 1 fm. (2.33) 

200 MeV ^ ' 

There is another point of view on the gluon condensate which I describe briefiy. In principle, 
all information about field theory is contained in the partition function being the functional 
integral over the fields. In the Euclidean formulation it is 

Z^JdA, exp (-^ / d'xF^^ ^-^ e-^^, (2.34) 
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where I have used that at large (Euchdean) time T the partition function picks up the ground 
state or vacuum energy S. For the sake of brevity I do not write the gauge fixing and Faddeev- 
Popov ghost terms. If the state is homogeneous, the energy can be written as £ = O^^V^^^ where 
6*^1, is the stress-energy tensor and V^^^ is the 3- volume of the system. Hence, at large 4- volumes 
V — V^^^T the partition function is Z — exp(— ^44!^). This ^44 includes zero-point oscillations 
and diverges badly. A more reasonable quantity is the partition function, normahzed to the 
partition function understood as a perturbative expansion about the zero-field vacuum^. 



exp 



-(^44-^44'-)^ ■ (2-35) 

We expect that the non-perturbative vacuum energy density ^44 — ^44''^' is a negative quantity 
since we have allowed for tunneling: as usual in quantum mechanics, it lowers the ground state 
energy. If the vacuum is isotropic, one has ^44 = 6'^^/4. Using the trace anomaly, 

where P{as) is the Gell-Mann-Low function, 

^ ^ _ _ _ (2.37) 



with 61^2 given by eq. (1.2), one gets [7]: 



■'P.T. 



f7^(i^.l/327r2)Npl (2.38) 



where (F^,^)np is the gluon field vacuum expectation value which is due to non-perturbative 
fluctuations, i.e. the gluon condensate. The aim of any QCD-vacTium builder is to minimize 
the vacuum energy or, equivalently, to maximize the gluon condensate. It is important that it 
is a renormalization-invariant quantity, meaning that its dependence on the ultraviolet cutoff 
II and the bare charge Q;s(/i) given at this cutoff is such that it is actually cutoff-independent. 
Such a combination is called A, see eq. (1.1). The gluon condensate has to be proportional to 
A^ by dimensions. 

The fact that the vacuum energy or, equivalently, the gluon condensate is a renormalization- 
invariant quantity leads to an infinite number of low-energy theorems [26] . Translated into the 
instanton- vacuum language, the renormalizability of the QCD implies that the probability that 
there are N Ps and /'s in the vacuum is [7, 27] 



P{N) ~ exp 




- 1 



(2.39) 



where (N) ~ y (F" )/(327r^) is the average number of J's and J's . 



^The latter can be distinguished from the former by imposing a condition that it does not contain integration 
over singular Yang-Mills potentials; recall that the instanton potentials are singular at the origins. 
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2.6 One-instanton weight 



The notion "instanton vacuum" implies that one assumes that the QCD partition function 
(2.34) is mainly saturated by an ensemble of interacting 7's and /'s , together with quantum 
fluctuations about them. Instantons are necessarily present in the QCD vacuum if only be- 
cause they lower the vacuum energy with respect to the purely perturbative (divergent) one. 
The question is whether they give the dominant contribution to the gluon condensate, and 
to other basic quantities. To answer this question one has to compute the partition function 
(2.34) assuming that it is mainly saturated by instantons, and to compare the obtained gluon 
condensate with the phenomenological one. 

The starting point of this calculation [7, 27] is the contribution of one isolated instanton 
to the partition function (2.34), or the one-instanton weight. We have already estimated the 
tunneling amplitude in eq. (2.18) but it is not sufficient: the prefactor is very important. To 
the 1-loop accuracy, it has been first computed by 't Hooft [3] for the SU{2) colour group, and 
generalized to arbitrary SU{N) by Bernard [28]. 

The general field can be decomposed as a sum of a classical field of an instanton A^(a;,^) 
where ^ is a set of 4A^c collective coordinates characterizing a given instanton (see eq. (2.28)), 
and of a presumably small quantum field a^{x): 

A^{x)^Aj^{x,0 + a,{x). (2.40) 

There is a subtlety in this decomposition due to the gauge freedom: an interested reader is 
addressed to ref. [7] where this subtlety is treated in detail. The action is 



Action = — j ct'xFl^ = ^ + — I cfx D^F^^a, + j d'^x a^W^^a, + O^o"). (2.41) 

Here the term linear in drops out because the instanton field satisfies the equation of motion. 
The quadratic form W^v has 4iVc zero modes related to the fact that the action does not depend 
on 4iVc collective coordinates. This brings in a divergence in the functional integral over the 
quantum field which, however, can and should be qualified as integrals over the collective 
coordinates: centre, size and orientations. Formally the functional integral over gives 



^det W^M')' 

which must be i) normalized (to the determinant of the free quadratic form, i.e. with no 
background field), ii) regularized (for example by using the Pauli-Villars method), and iii) 
accounted for the zero modes. Actually one has to compute a "quadrupole" combination, 

' dei'W det(iyo + /i^) " 
detWo det(iy + //2) 

where Wq is the quadratic form with no background field and //^ is the Pauli-Villars mass 
playing the role of the ultraviolet cutoff; the prime reminds that the zero modes should be 



(2.43) 
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removed and treated separately. The resulting one-instanton contribution to the partition 
function (normalized to the free one) is [3, 28]: 



■inst 



-2^P.T. 

(io(p) 



27T 



2Nc 



as{p) 



exp 



27r 



«s(/^). 



(2.44) 
(2.45) 



The fact that there are all in all 47Vc integrations over the collective coordinates z^, p, C/ re- 
flects ANc zero modes in the instanton background. The numerical coefficient C{N^ depends 
implicitly on the regularization scheme used. In the Pauli-Villars scheme exploited above [28] 



C(7Vc) 



4.60exp(-1.68A^c) 



7r2(7V^-l)!(Ar^-2)!- 
If the scheme is changed, one has to change the coefficient C{Nc) C'{Nc 
One has [29]: Ap.y. = e^K^^ = 40.66 e A 



(2.46) 

C{N,) ■ (A/A')^ 



lat 



Eq. (2.45) cannot yet be expressed through the 2-loop renormalization-invariant combination 
A (1.1) as it is written to the 1-loop accuracy only. In the 2-loop approximation the instanton 
weight is given by [30, 7] 



doip) = ^Pipr^ 



exp 



-P'\p)+i2N, 



b2 \ b2 ln/3(p) 
'2bj2bi (3{p) 



O 



Mp), 



(2.47) 



■(AP)- 



where (3{p) = 2n/as{p) and /?^^(p) are the inverse charges to the 1-loop and 2-loop accuracy, 
respectively (not to be confused with the Gell-Mann-Low function!): 



(2.48) 

34 

b2 = yA^e'- (2-49) 

These equations express the one-instanton weight do{p) through the cutoff-independent combi- 
nation A (1.1), and the instanton size p. This is how the 'transmutation of dimensions' occurs 
in the instanton calculus and how A enters into the game. Henceforth all dimensional quantities 
will be expressed through A, which is very much welcome. 

Notice that the integral over the instanton sizes in eq. (2.44) diverges as a high power of p 
at large p: this is of course the consequence of asymptotic freedom. It means that individual 
instantons tend to swell. This circumstance plagued the instanton calculus for many years. If 
one attemts to cut the p integrals "by hand" , one violates the renormalization properties of the 
YM theory, as mentioned in the previous section. Actually the size integrals appear to be cut 
from above due to instanton interactions. 



m 



62 , 2/3 p 



bi In 
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2.7 Instanton ensemble 



To get a volume effect from instantons one needs to consider an II ensemble, with their total 
number N proportional to the 4-dimensional volume V . Immediately a mathematical difficulty 
arises: any superposition of /'s and J's is not, strictly speaking, a solution of the equation of 
motion, therefore, one cannot directly use the semiclassical approach of the previous section. 
One way to overcome this difficulty is to use a variational principle [7]. Its idea is to use a 
modified YM action for which a chosen II ansatz is a saddle point. Exploiting the convexity 
of the exponent one can prove that the true vacuum energy is less than that obtained from 
the modified action. One can therefore use variational parameters (or even functions) to get a 
best upper bound for the vacuum energy. It is not the Rayleigh-Ritz but rather the Feynman 
variational principle since the method has been suggested by Feynman in his famous study 
of the polaron problem. The gauge theory is more difficult, though: one has not to loose 
either gauge invariance or the renormahzation properties of the YM theory. These difficulties 
were overcome in ref. [7], see also [27]. It should be kept in mind that we are dealing with 
"strong interactions", meaning that all dimenionless quantities are generally speaking of the 
order of unity - there are no small parameters in the theory. Therefore, one has to use certain 
approximate methods, and the variational principle is among the best. Todays direct lattice 
investigation of the // ensemble seem to indicate that we have obtained rather accurate numbers 
in this difficult problem. 

In the variational approach, the normalized (to perturbative) and regularized YM partition 
function takes the form of a partition function for a grand canonical ensemble of interacting 
pseudoparticles of two kind, /'s and /'s : 

Z 1 1 ^++^- f 

H WlWl n d'^ZndpndUi do{pn) exp{-Uint), (2.50) 

where (io(p) is the 1-instanton weight (2.47). The integrals are over the collective coordinates of 
(anti)instantons: their coordinates 2;, sizes p and orientations given by SU{Nc) unitary matrices 
U; dU means the Haar measure normalized to unity. The instanton interaction potential f/jnt 
(to be discussed below) depends on the separation between pseudoparticles, Zm — Zn, their sizes 
Pm,n and their relative orientations UmUl- In the variational approach the interaction between 
instantons arise from i) the defect of the classical action, ii) the non-factorization of quantum 
determinants and Hi) the non-factorization of Jacobians when one passes to integration over the 
collective coordinates. All three factors are ansatz-dcpcndent, but there is a tendency towards 
a cancellation of the ansatz-dependent pieces. Qualitatively, in any ansatz the interactions 
between 7's and Ps resemble those of molecules: at large separations there is an attraction, at 
smaller separations there is a repulsion. It is very important that the interactions depend on 
the relative orientations of instantons: if one averages over orientations (which is the natural 
thing to do if the // medium is in a disordered phase; if not, one would expect a spontaneous 
breaking of both Lorentz and colour symmetries [7]), the interactions seem to be repulsive at 
any separations. 

In general, the mere notion of the instanton interactions is notorious for being ill-defined 
since instanton -|- antiinstanton is not a solution of the equation of motion. Such a configuration 
belongs to a sector with topological charge zero, thus it seems to be impossible to distinguish 
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it from what is encountered in perturbation theory. The variational approach uses brute force 
in deahng with the problem, and the results appear to be somewhat dependent on the ansatz 
used. Thanks to the inequality for the vacuum energy mentioned above, we still get quite a 
useful information. However, recently a mathematically unequivocal definition of the instanton 
interaction has been suggested, based on the one hand on analyticity and unitarity [31] and on 
the other hand on certain singular solutions of the YM equations of motion [32]. Both definitions 
cut off automatically contributions of the perturbation theory. The first three leading terms 
for the interaction potential at large separations has been computed by the two very different 
methods [31, 32] with coinciding results. At smaller separations one observes a strong repulsion 
[32]. 

At this point I should mention certain experience one gains from a simpler 2-dimensional so- 
called CP^ model, also possessing instantons as classical Euclidean solutions. Contrary to the 
4d YM theory, the instanton measure in that model is known exactly [33, 34]. In the dilute limit 
the instanton measure reduces to the product of integrals over instanton sizes, positions and 
orientations, as in eq. (2.50). The exact measure, however, is written in terms of the so-called 
'instanton quarks' which does not suppose that instantons are dilute. The statistical mechanics 
of J's and J's in this model has been studied in ref. [35] both by analytical methods and by 
numerical simulations. Although the 'instanton quark' parametrization allows for complete 
'melting' of instantons and is quite opposite in spirit to the dilute-gas ansatz, it has been 
observed that, owing to a combination of purely geometric and dynamic reasons, the vast 
majority of 'instanton quarks' form neutral clusters which can be identified with well-separated 
instantons. Of course, there is always a fraction of overlapping instantons in the vacuum, 
however, it is small even in the 2d case; in the Ad YM case both reasons mentioned above are 
expected to be even stronger. 

Summing up the discussion, I would say that today there exists no evidence that a variational 
calculation with the simplest sum ansatz used in ref. [7] is qualitatively or even quantitatively 
incorrect, therefore I will cite the numerics from those calculations in what follows. The main 
finding [7, 27] is that the // ensemble (2.50) stabilizes at a certain density related to the A 
parameter (there is no other dimensional quantity in the theory!) 



The average instanton size and the average separation between instantons are, respectively. 



(F^2,/32vr2) ^ -(Q^) > (0.75A^)^ (2.51) 



p ~ 0.48/ A^jg ~ 0.35 fm, (2.52) 
R = i^yj ~1.35/Aj^~0.95fm, (2.53) 

if one uses Aj^ = 280 MeV as it follows from the DIS data. Earlier, very similar charactersitics, 
P — I fm, ^ = 1 fm, have been suggested by Shuryak [5] from studying the phenomenological 
applications of instantons. 

Instanton interactions lead to the modification of the (divergent) size distribution function 
do{p) (2.47) by a distribution decreasing at large p. The use of the variational principle yields 
a Gaussian cutoff for large sizes [7, 27]: 
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Figure 3: "Cooling" the normal zero-point oscillations reveals large fluctuations of the gluon 
field, which were identified with instantons and anti-instantons with random positions and 
sizes [21]. The left column shows the action density and the right column shows the topological 
charge density for the same snapshot. 



In fact, it is a rather narrow distribution peaked around p (2.52); therefore for practical esti- 
mates in what follows I shall just replace all instantons by the average-size one. 

It should be said that, strictly speaking, nothing can prevent some instantons to be anoma- 
lously large and overlapping with other. For overlapping instantons the notion of size distribu- 
tion becomes senseless. The question is quantitative: how often and how strong do instanstons 
overlap. Given the estimate (2.52,2.53), it seems that the majority of instantons in the vacuum 
ensemble are well-isolated. 

In the recent years instantons have been intensively studied by direct numerical simulations 
of gluon fields on the lattice, using various configuration-smoothing methods [21, 36, 37]. A 
typical snapshot of gluon fluctuations in the vacuum is shown in Fig. 3 borrowed from ref. [21]. 
Naturally, it is heavily dominated by normal perturbative UV-divergent zero-point oscillations 




(2.54) 
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of the field. However, after "cooling" down these oscillations one reveals a smooth background 
field which was shown in rcf. [21] to be nothing but an ensemble of instantons and anti- 
instantons with random positions and sizes The lower part of Fig. 3 is what is left of the 
upper part after "cooling" that particular configuration. The average sizes and separations of 
instantons found vary somewhat depending on the concrete smearing method used. Ref. [21] 
gives the following values 

p ~ 0.36 fm, R = {N/V)--i ~ 0.89 fm, (2.55) 
which are not far from the estimate from the variational principle. The ratio, 

seems to be more stable: it follows from phenomcnological [5], variational [7, 27] and lattice 
[21, 36, 37] studies. It means that the packing fraction, i.e. the fraction of the 4-dimensional 
volume occupied by instantons appears to be rather small, vr^p*^/^^ ~ 1/8. This small packing 
fraction of the instantons gives an a posteriori justification for the use of the semi-classical 
methods. As I shall show in the next sections, it also enables one to identify adequate degrees 
of freedom to describe the low-energy QCD. 

3 Non-instanton semiclassical configurations 

Instantons induce certain potential between static quarks in the pure glue theory, which de- 
pends on the instanton size distribution [4, 14, 38]. For fast convergent size distributions like 
that given by eq. (2.54) the potential first rises as function of the interquark separation but 
asymptotically it flattens out. If the size distribution happens to fall off as d{p) ~ 1/p^ at large 
p one gets a linear infinitely rising potential [38]. However, such a size distribution means that 
large instantons inevitably overlap, which means that the "center, size, orientations" collective 
cooredinates are not the most adequate. Standard instantons do not induce an infinitely rising 
linear potential between static quarks in the pure glue theory. This observation stimulates 
the search for other semiclassical gluonic objects that could be responsible for confinement. 
Among them are merons [39, 40], calorons with non-trivial holonomy [41, 42], BPS monopoles 
or dyons [43, 44], vortices [45, 46], etc. I describe briefiy these objects below. 

Merons 



Meron is a self-dual solution of the YM equation of motion [39] , 

^:^=%^'^^ (3-1) 

which is interesting because at fixed time its spatial components reminds the field of a 3d 
monopole of size t: 

^Quite recently a more involved fluctuation-smearing procedure carried on the lattice has indicated that 
instantons might have an additional structure, see the next section. 
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A'^ — F ■ ■ 



(3.2) 



However, its action diverges logarithmically both in the UV and IR regions; therefore such 
configurations are not encountered in the vacuum. Callan, Dashen and Gross [40] suggested to 
consider a regularized meron pair, 



Aa 



V. 



1X1/ 



[X 



+ 



[X 



Z2 



(3.3) 



{x-ziY + pl {x 

whose action is finite, and at \x\ — > oo it becomes the field of an instanton. Eq. (3.3) corresponds 
to the field of a monopole pair which has been created, moved to a separation R — \zi — Z'z] 
and then annihilated after the time R. Merons can be viewed as "half-instantons" . The idea 
of Callan, Dashen and Gross was that instantons could dissociate into a gas of (regularized) 
merons, which could lead to confinement. Recently, a relation of merons to monopoles and 
vortices has been discussed - see ref. [47] where such relation has been studied numerically, 
and references to earlier work therein. So far merons have not been identified from lattice 
configurations. Although the idea that meron configurations may well prove to be relevant to 
confinement in 4(i, I have certain reservations as confinement is also a property of the 3(i pure 
YM theory where there are no merons. Therefore, if confinement is due to some semiclassical 
"objects" , 3d and 2d field configurations may be of more relevance. 



Dyons and calorons with non-trivial holonomy 



Dyons or Bogomolnyi-Prasad-Sommerfeld (BPS) monopoles are self-dual solutions of the YM 
equations of motion with static {i.e. time-independent) action density, which have both the 
magnetic and electric field at infinity decaying as Therefore these objects carry both 

electric and magnetic charges. In the 3 -|- 1- dimensional SU{2) gauge theory there are in fact 
two types of self-dual dyons [48]: M and L with (electric, magnetic) charges (-I-, -|-) and (— , — ), 
and two types of anti-self-dual dyons M and L with charges (-I-, — ) and (— , -|-), respectively. 
Their explicit fields can be found e.g. in ref. [49]. In SU{Nc) theory there are 2Nc different 
dyons [48, 50]: Mi, M2, ...Mjv^_i ones with charges counted with respect to A^c — 1 Gartan 
generators and one L dyon with charges compensating those of Mi...M^^^i to zero, and their 
anti-self-dual counterparts. 

Speaking of dyons one implies that the Euclidean space-time is compactified in the 'time' 
direction whose inverse circumference is temperature T, with the usual periodic boundary 
conditions for boson fields. However, the temperature may go to zero, in which case the 4d 
Euclidean invariance is restored. 

Dyons' essence is that the component of the dyon field tends to a constant value at 
spatial infinity. This constant /I4 can be eliminated by a time-dependent gauge transformation. 
However then the fields violate the periodic boundary conditions, unless A^ has quantized 
values corresponding to trivial holonomy, see below. Therefore, in a general case one implies 
that dyons have a non-zero value of at spatial infinity. 

A single dyon can be considered in whatever gauge. The simplest is the "hedgehog" gauge 
in which the dyon field is smooth everywhere. However, if one wishes to have a vacuum filled 
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by dyons, one has to take more than one dyon. Two and more dyons can be put together 
only if all colour components of at infinity are the same for all dyons. Before assembling 
"hedgehogs" together, their "hair" has to be "gauge-combed" in such a way that A^t" at infinity 
is a constant matrix, which can well be chosen to be diagonal, i.e. belonging to the Cartan 
subalgebra. The gauge meeting this requirement is called "stringy" gauge as it necessarily has 
a pure-gauge string-like singularity starting from the dyon center. For explicit formulae see e.g. 
the Appendix of rcf. [49]. In other words two or more dyons need to have the same orientation 
in colour space. This orientation is preserved throughout the 3d volume. The A4 component 
of the YM field plays here the role of the Higgs field in the adjoint representation. The spatial 
size of M dyons is l/\Al\ (for the SU{2) group) and that of L dyons is l/|27rr - Al\ where 
T is the Euclidean temperature. The mere notion of the ensemble of dyons (or monopoles) 
implies that colour symmetry is in a sense spontaneously broken. Of course, once colour is 
aligned, one can always randomize the colour orientation by an arbitrary point- dependent 
gauge transformation, just as the direction of the Higgs field can be randomized but that does 
not undermine the essence of the Higgs effect. 

This is very distinct from instantons which can have arbitrary colour orientations in the 
ensemble. Mathematically, one can see it from the number of zero modes. For example, for 
the SU (2) group instantons have 4 translation, 1 dilatation and 3 orientation zero modes, see 
subsection 2.4. A dyon (in any gauge group) has only 4 zero modes out of which 3 are spatial 
translations and 1 has a double description: one can call it the spurious translation in the 'time' 
direction or, better, a global U{1) phase. The dyon solution has no collective coordinates that 
can be identified with colour orientation. This fact is not accidental but related precisely to 
that dyons have a non- vanishing field at infinity. One can formally introduce an 'orientation' 
degree of freedom but the corresponding zero mode will be not normalizable. In the instanton 
case all fields decay fast enough at infinity to make the orientation mode normalizable and 
hence physical. 

The fact that dyons have non-zero A4 at spatial infinity means also that the Polyakov line 
(also called the holonomy) is in general non-trivial at infinity, i.e. does not belong to the group 
center, unless assumes certain quantized values: 



Evaluating the effects of quantum fiuctuations about a dyon proved to be a difficult task, 
much more difficult than that about instantons. The most advanced calculation is in ref. [51] 
where quantum determinants have been computed in some limiting case. Probably the tech- 
nique developed therein allows one to read off the result for determinants in a general case but 
this has not been done. One thing is clear, however, without calculations: the quantum deter- 
minants about dyons strongly diverge in the IR region! The point is, at non-zero temperatures 
the quantum action contains, inter alia., a potential-energy term [52, 53] 




(3.4) 




[for ,SC/(2) group]. 



(3.5) 



As follows from eq. (3.4) the trace of the Polyakov line is related to (j) as 
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iTrP^cos;^. (3.6) 



2T 

The zero energy of the potential corresponds to P = ±1, i.e. to the trivial holonomy. If a 
dyon has ^ 2TrTn at spatial infinity the corresponding quantum action is positive-definite 
and proportional to the 3d volume. Therefore, dyons with non-trivial holonomy seem to be 
strictly forbidden: quantum fluctuations about them have an unacceptably large action! 



There are two known generalizations of instantons at non-zero temperatures. One is the 
periodic instanton of Harrington and Shepard [54] studied in detail in ref. [52]. These periodic 
instantons, also called calorons, have trivial holonomy at spatial infinity, therefore they are not 
suppressed by the above quantum potential energy. The vacuum made of those instantons has 
been investigated, using the variational principle, in ref. [55]. 

The other generalization has been constructed recently by Kraan and van Baal [41] and 
Lee and Lu [42]; it has been named caloron with non-trivial holonomy. I shall call it for short 
the KvBLL caloron. It is also a self-dual solution of the YM cqTiations of motion with a unit 
topological charge. The fascinating feature of this construction is that it can be viewed as 
"made of" one L and one M dyon, with total zero electric and magnetic charges. Although 
the action density of isolated L and M dyons does not depend on time, their combination in 
the KvBLL solution is generally non-static: the L, M "consituents" show up not as 3d but 
rather as 4d lumps. When the temperature goes to zero, these lumps merge, and the KvBLL 
caloron reduces to the usual instanton (as does the standard Harrington-Shepard caloron), plus 
corrections of the order of T. However, the holonomy remains fixed and non-trivial at spatial 
infinity. It means that quantum fiuctuations strongly suppress individual KvBLL calorons at 
any temperature, as they suppress single L, M dyons. 

It is very interesting that precisely these objects determine the physics of the supersymmetric 
YM theory where in addition to gluons there are gluinos, i.e. Majorana (or Weyl) fermions 
in the adjoint representation. Because of supersymmetry, the boson and fermion determinants 
about L, M dyons cancel exactly, so that the perturbative potential energy (3.5) is identically 
zero for all temperatures (actually to all loops). Therefore, in the supersymmetric theory dyons 
are openly allowed. [To be more precise, the cancellation occurs when periodic conditions for 
gluinos are imposed, so it is the compactification in one (time) direction that is implied, rather 
than physical temperature which requires antiperiodic fermions. However, I shall still call 
it "temperature".] Further on, it turns out [56, 50] that dyons generate a non-perturbative 
potential having a minimum ai (f) — ttT, i.e. where the perturbative potential would have the 
maximum. This value of A4 corresponds to the holonomy Tr P = at spatial infinity, which is 
the "most non-trivial"; as a matter of fact it is one of the confinement's requirements. However, 
it implies that A4 has to lie in some direction in colour space thus breaking the colour group 
to the maximal Abelian subgroup, at least at small compactification circumference [49]. 

In the supersymmetric YM theory there is a non-zero gluino condensate. It is analogous 
to the quark condensate in QCD, which we consider in the next section. However, contrary to 
QCD, in the supersymmetric case the gluino condensate can be computed exactly and expressed 
- via the transmutation of dimensions - through the scale A. Its correct value is reproduced 
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by saturating the gluino condensate by L, M dyons' zero fermion modes [56]. On the contrary, 
the saturation of the (square of) gluino condensate by instanton zero modes gives the wrong 

result, namely y^4/5 that of the correct value [57]. The paradox is that both derivations are 
seemingly clean. 

This striking 20-year-old paradox has been recently resolved [49] by the observation that 
the (square of) gluino condensate must be computed not in the instanton background but in 
the background of exact solutions "made of" LL, MM and LM dyons. The first two are the 
double-monopole solutions and the last one is the KvBLL caloron. As the temperature goes to 
zero, the LL and MM solutions have locally vanishing fields, whereas the KvBLL LM solution 
reduces to the instanton field up to a locally vanishing difference. Therefore, naively one would 
conclude that the dyon calculation of the gluino condensate, which is a local quantity, should 
be equivalent to the instanton one, but it is not! The fields vanishing as the inverse size of the 
system have a finite effect on such a local quantity as the gluino condensate! This is quite un- 
usual. The crucial difference between the (wrong) instanton and the (correct) dyon calculations 
is in the value of the Polyakov loop, which remains finite. In the supersymmetric YM theory 
configurations having Tr P = at infinity are not only allowed but dynamically preferred as 
compared to those with P = ±1. In non-supersymmetric theory it looks as if it was the opposite. 

An intriguing question is whether the pcrturbative potential energy (3.5) which forbids 
individual dyons in the pure YM theory can be overruled by non-perturbativc contributions of 
an ensemble of dyons. I think that there is such a chance. Let us consider the case of non- 
zero temperatures, and let us imagine an ensemble of L, M, Z, M dyons. Assuming that the 
temperature is high enough one can suppose that dyons are dilute and hence i) it is actually 
an ensemble of static objects, ii) their interaction can be neglected as compared to their own 
action. The M,M action is whereas the L, L action is ^■'^^^^ where = ^/AfAl 

and Al is the common value of the dyons' field at spatial infinity. The weight of one dyon in 
the semiclassical approximation is exp(— Action), times the integral over collective coordinates 
/ d^z, times something which makes it dimensionless - usually it is the solution's size. Therefore, 
I would write the partition function of the dyon gas as 
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(3.7) 



This is naturally a very crude estimate. Interactions have been neglected. Exact quantum 
determinants are unknown. However, it is known that they make the gauge coupling run, i.e. 
their role is basically to replace the coupling constants in eq. (3.7) by the running coupling at the 

temperature scale: ^ Ittt j ' "^^^ overall numerical factor cannot be determined from 

such simple considerations - it must be computed. When evaluating quantum determinants 
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Figure 4: An illustration of how dyons can induce the confinement-deconfinement phase transi- 
tion. We plot the sum of the perturbative (3.5) and dyon-induced (3.8) potentials as function of 
0/27rT. We choose the constants A = 280 MeV, c = 2. The curves corresponds to temperatures 
T = 250 MeV (sohd), T = 280 MeV (long-dashed) and T = 400 MeV (short-dashed). In this 
example, the phase transition is at Tc = 280 MeV. 



one implies that they are normalized not to the free ones but those in constant background 
field of A4, to make the normalized determinants IR finite. The determinants in the constant 
field are known to produce the perturbative potential (3.5). 

An educated guess for the non-perturbative potential induced by dyons is therefore 
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to which the perturbative potential (3.5) must to be added. Let me note that in the supersym- 
metric theory there is a similar non-perturbative potential [56] (but in that case it is an exact 
result) and the perturbative potential is absent. 

The perturbative potential (3.5) has minima at = 0, 27rT corresponding to trivial holon- 
omy P = ±1, while the non-perturbative potential (3.8) has the minimum just in the middle, at 
= ttT corresponding to the non-trivial holonomy Tr P = 0, which is one of the requirements 
of confinement. Clearly, at very large temperatures the perturbative potential wins, and the 
system settles in the perturbative vacuum with trivial holonomy. At temperatures below certain 
critical Tc = const. A the non-perturbative potential prevails, dyons with the asymptotic value 
of = 7tT are preferred, and the system presumably goes into the confinement phase. For the 
particular model of the dyon-induced potential (3.8) it is a second-order phase transition, see 
Fig. 4. For SU{3) and higher groups both the perturbative and non-perturbative potentials 
depend on more variables than just 0, and one can well imagine their interplay leading to a 
first order transition. I have presented here what might be a crude microscopic model of the 
confinement-deconfinement transition, based on dyons. 

It should be noted that recent lattice simulations support the presence of KvBLL calorons 
(or, more generally, dyons) near the phase transition point [58, 59]. 
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As the temperature lowers, the ensemble of L, L, M, M dyons becomes more dense, and 
these objects (which have time-independent action densities when they are well-isolated) do 
not represent a static field anymore: the action density will have 4d lumps. Their statistical 
mechanics must be most intriguing! As I mentioned before, at T — > the KvBLL calorons 
tend to instantons, with their L, M 'constituents' merged together. This is the most dehcate 
moment. I do not know what prevents them from glueing up completely into normal instantons, 
and what guarantees that non-trivial holonomy is dynamically preferred. 

There are problems with the dyon scenario at high temperatures as well. Dyons survive 
for a while above Tc although the minima of the potential shown in Fig. 4 correspond to the 
holonomy approaching P — ±1. Eventually, the Ad theory at high temperatures becomes the 
3d theory. We beheve that confinement is also a property of the pure-glue 3d YM theory but 
in the true 3d case there are no dyons as there is no A4 component of the field, which is critical 
for the dyon construction, and we do not know of any other classical solutions there. 

Monopoles and vortices 

A monopole is, by definition, a 3d field configuration whose magnetic field is decaying as 
at large distances. A center vortex is a 2d field configuration such that a Wilson loop winding 
about it takes the value of one of the non-trivial elements of the group center Z(Ac)- A simple 
dimensional analysis based on a rescaling of the assumed solutions' size shows that there are 
no such solutions of the classical YM equations apart from BPS dyons described above but in 
the true 3d case they are absent. It should be noted however that many non-trivial classical 
solutions appear when one imposes twisted instead of periodic boundary conditions [60]. 

In the lattice community a completely different definition of monopoles and vortices is used, 
for practical reasons. They are defined through gauge fixing, see a recent review by Greensite 
[61]. Monopoles and vortices identified on the lattice do not necessarily correspond to any 
semiclassical field configurations. Lattice monopoles are points in 3d or lines in Ad space while 
lattice center vortices are lines in 3d and surfaces in Ad. If there is any materialistic meaning 
behind them both monopoles and vortices need to be 'thick', that is to have on the average a 
finite radius of the order of 1/A since it is the only physical scale in the YM theory. This hints 
that they need not be classical solutions (the more so that there are none) but rather local 
minima of the effective action with quantum corrections included. 

To give an example, in high-temperature 2+1 dimensional pure YM theory a 'thick' center 
vortex solution has beed found as a local minimum of the classical plus 1-loop quantum ac- 
tion [62] . In 3d and Ad the 1-loop effective action about a vortex has been computed in ref . [63] 
and the corresponding local minima found, too. So far no similar analysis has been performed 
for 3d monopoles, which would be of some interest. Although the radius of the vortex has been 
found in units of 1/A, the size appears to be too large and the minimum is very shallow. Later 
on it was discovered that vortices have a negative mode [64] i.e. are unstable. These studies are 
rather discouraging for attempts to attribute any semiclassical meaning to vortices. It should 
be added that nothing is known about the statistical mechanics of vortices from the theory 
side, nor even is the object itself defined. 

The idea that some semiclassical objects may be relevant for confinement faces certain 
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difficulty of a general nature. It is widely believed that the string tension between static quarks 
is stable in the limit of large number of colours Nc- This is supported by lattice observations 
at medium Nc [65, 66]. Assuming the simplest scenario that confinement is driven by weakly 
interacting center vortices populating the YM vacuum, one finds then that the density of 
vortices piercing any given 2d plane should rise linearly with Nc [67]. It means that at large 
enough Nc vortices have to overlap. What happens then is not clear. If they remain independent 
despite overlap, it should be clarified how it happens. If they are not independent but interact 
strongly, the simple argument why vortices lead to the area behaviour of large Wilson loops is 
lost. 

In fact all semiclassical objects face a problem at large Nc, and the reason is quite general: 
they possess a finite action density Tr F^^. If semiclassical objects saturate the non-perturbative 
gluon condensate (see subsection 2.5) their density has to rise linearly with N^, since the gluon 
condensate does so, as it follows from Nc counting rules. This is a universal problem for 
instantons, dyons and vortices: they have to overlap at large Nc- A way to prevent these 
objects from complete "melting" is as follows. Start with instantons. The use of the variational 
principle [7, 27] confirms the expectation that instanton density is 0{Nc) whereas their average 
size is 0(1), i.e they have to overlap at large Nc- However, instanton is basically an SU{2) 
object, and it has been noticed by T. Schafer [68] that instantons residing in different "corners" 
of the large SU {N) group space do not "talk" to each other, meaning that they might well 
be statistically independent. A more precise meaning of it is suggested by dyons. Same as 
instantons, dyons have finite action and finite size, however there are precisely Nc different 
types of dyons (see above), therefore each type of dyons has a finite density at large Nc although 
different types have to overlap. A dyon is an SU{2) object built about one of the Cartan 
generators of the SU{N) group, diag(0...0, 1, — 1, 0...0). Clearly, only dyons that are nearest 
neighbours in colour space do interact, however 0{Nc) types of dyons are transparent for any 
given type. Therefore, there is no harm if dyons of the not-nearest-neighbour types sit on top 
of each other. Dyons from commuting SU{2) subgroups do not interact either in the classical 
or in the quantum sense. Probably a similar argument can be given for vortices, as there are 
Nc types of them. 

Thus, the transition to infinite Nc in terms of semiclassical objects may be in fact quite 
smooth. It is very important to understand this transition clearly since lattice simulations 
indicate [65, 66] that already Nc — 2 is not so far away from infinity. Apparently, confinement 
should be explained by the same mechanism both at small and large Nc- Also, the asymptotic 
area law for large spatial Wilson loops is believed to be the property of all theories starting 
from 2-1-1 dimensions and to 4 dimensions, and it would be natural if it is driven basically by 
the same mechanism. 

At the same time it should be stressed again that at large Nc any semiclassical picture 
implies that objects are very dense, and hence their colour fields are fast varying. It indicates 
that the semiclassical methodology is probably not the best way to understand confinement. 
Personally, I think that passing to dual gauge-invariant variables [69] which do not need to 
oscillate violently, is a more promising strategy in solving this awful problem. 
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4 Chiral symmetry breaking by instantons 



4.1 Chiral symmetry breaking by definition 

The QCD Lagrangian with Nf massless flavours is known to posses a large global symmetry, 
namely a symmetry under U{Nf) x U{Nf) independent rotations of left- and right-handed 
quark fields. This symmetry is called chiral Instead of rotating separately the 2-component 
Weyl spinors corresponding to left- and right-handed components of quark fields, one can make 
independent vector and axial U (Nf) rotations of the full 4-component Dirac spinors - the QCD 
lagrangian is invariant under these transformations too. 

Meanwhile, axial transformations mix states with different P-parities. Therefore, were that 
symmetry exact, one would observe parity degeneracy of all states with otherwise the same 
quantum numbers. In reality the splittings between states with the same quantum numbers 
but opposite parities are huge. For example, the splitting between the vector p and the axial 
Oi meson is (1260 — 770) ~ 500 MeV; the splitting between the nucleon and its parity partner 
is even larger: (1535 - 940) ~ 600 MeV. 

The splittings are too large to be explained by the small bare or current quark masses which 
break the chiral symmetry from the beginning. Indeed, the current masses of light quarks are: 
m„ ~ 4 MeV, rrid ~ 7 MeV, ~ 150 MeV. The conclusion one can draw from these numbers 
is that the chiral symmetry of the QCD Lagrangian is broken down spontaneously, and very 
strongly. Consequently, one should have light (pseudo) Goldstone pseudoscalar hadrons - their 
role is played by pions which indeed are by far the lightest hadrons. 

The order parameter associated with chiral symmetry breaking is the so-called chiral or 
quark condensate: 



It should be noted that this quantity is well defined only for massless quarks, otherwise it is 
somewhat ambiguous. By definition, this is the quark Green function taken at one point; in 
momentum space it is a closed quark loop: 



If the quark propagator is massless and has only the 'slash' term, the trace over the spinor 
indices in the loop gives an identical zero. Therefore, chiral symmetry breaking implies that a 
massless (or nearly massless) quark develops a non-zero dynamical mass M{k), i.e. a 'non-slash' 
term in the propagator. There are no reasons for this quantity to be a constant independent 
of the momentum; moreover, we understand that it should anyhow vanish at large momentum. 
Sometimes it is called the constituent quark mass, however a momentum-dependent dynamical 
quark mass M{k) is a more adequate term which I shall use below. 

The spontaneous generation of the dynamical quark mass (equivalent to the spontaneous 
chiral symmetry breaking, SCSB) is the most important feature of QCD being key to the whole 
hadron phenomenology. The theory's task is to get M{k) in the form 

^The word was coined by Lord Kelvin in 1894 to describe moleculas not superimposable on its mirror image. 
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(4.2) 



26 



M{k) = Af{k/A) (4.3) 

where A is the renormaUzation-invariant combination (1.1) and / is some function. Instantons 
enable one to get M{k) in the needed form and to find the function. But first let us derive 
some general relations. 

We start by writing down the QCD partition function. Functional integrals are well defined 
in Euclidean space which is obtained by the following formal substitutions of Minkowski space 
quantitites: 



IX MO 



XE4:, 

4, 



XMi 

7mo 



lEA, 



Amo 

iMi = 



- i^EA, A Mi — Aei, 



(4.4) 



Neglecting for brevity the gauge fixing and Faddeev-Popov ghost terms, the QCD partition 
function with quarks can be written as 



= I DA^Dil^Dil^^ exp j + j ^^Mf 

= / DA^ exp 



/ 



ndet(iy+im/). 

/ 



(4.5) 



The chiral condensate of a given fiavour / is, by definition. 

Id 



The Dirac operator has the form 



if^^^{id^ + A'^ + ai,) 



(4.6) 



(4.7) 



where A^^ denotes the classical field of the // ensemble and is a presumably small field of 
quantum fluctuations about that ensemble, which I shall neglect as it has little impact on chiral 
symmetry breaking. Integrating over DA^^ in eq. (4.5) means averaging over the // ensemble 
with the partition function (2.50), therefore one can write 



det(iy + im) 



(4.8) 



where I temporarily restrict the discussion to the case of only one fiavour for simplicity. Because 
of the im term the Dirac operator in (4.8) is formally not Hermitian; however the determinant is 
real due to the following observation. Suppose we have found the eigenvalues and eigenfunctions 
of the Dirac operator. 



(4.9) 
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then for any An 7^ there is an eigcnfunction $„/ = 75$^ whose eigenvalue is A^/ = — A„. This 
is because 75 anticommutes with iV. Owing to this the fermion determinant can be written as 



det(iy + im) 



n(An + im) = \ll[{Xl + m 

n 



exp 



exp 



- dXiy{X)\n{X'' + m^) , z/(A) = ^ 5(A - A,), (4.10) 



where I have introduced the spectral density z/(A) of the Dirac operator zV. Note that the 
last expression is real and even in m, which is a manifestation of the QCD chiral invariance. 
Differentiating eq. (4.10) in m and putting it to zero one gets according to the general eq. (4.6) 
a formula for the chiral condensate: 



^JdX iy{X) ln(A2 + m") 



V dm 

1 f°° m 



1 fOO 



A^ + 



m— ►0 

(4.11) 



where z/(A) means averaging over the instanton ensemble together with the weight given by the 
fermion determinant itself. The latter, however, may be cancelled in the so-called quenched 
approximation where the back influence of quarks on the dynamics is neglected. Theoretically, 
this is justified at large Nc- Naively, one would think that the r.h.s. of eq. (4.11) is zero at 
m — s> 0. That would be correct for a finite-volume system with a discrete spectral density. 
However, if the volume goes to infinity faster than m goes to zero (which is what one should 
assume in the thermodynamic limit) the second factor in the integrand becomes a representation 
of a 5-function, 



'"^->°sign(m)7r5(A), (4.12) 



A^ + m^ 



{jjjjjj) = -— sign(m)7ri/(0). (4.13) 



so that one obtains [10] : 



It is known as the Banks-Casher relation [70]. The chiral condensate is thus proportional to 
the averaged spectral density of the Dirac operator at zero eigenvalues. 

The appearance of the sign function is not accidental: it means that at small m QCD 
partition function depends on m non-analytically: 



\nZ = V{co + 7rz/(0)|m| + C2m^\n{\m\) + ...). (4.14) 

The fact that the partition function is even in m is the refiection of the original invariance of the 
QCD under 75 rotations; the fact that it is non-analytic in the symmetry-breaking parameter 
m is typical for spontaneous symmetry breaking. A generalization of the above formulae to the 
case of several quark flavours can be found in ref . [71] . 
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4.2 Physics: quarks hopping from one instanton to another 



Below I follow refs. [10, 11]. 

The key observation is that the Dirac operator in the background field of one (anti) instanton 
has an exact zero mode with A = [3]. It is a consequence of a general Atiah-Singer index 
theorem; in our case it is guaranteed by the unit Pontryagin index or the topological charge 
of the instanton field. These zero modes are 2-component Weyl spinors: n^'/i^-handed for 
instantons and /e/f-handed for antiinstantons. Explicitly, the zero modes are {a — l...Nc is the 
colour and i,j,k — 1, 2 are the spinor indices): 

(f){x, p) 

Here 2:1^, pi, Ui are the center, size and orientation of an instanton and Z21J,, P2, U2 are those of 
an antiinstanton, respectively, e-'^ is the 2x2 antisymmetric matrix. 

For infinitely separated / and / one has thus two degenerate states with exactly zero eigen- 
values. As usual in quantum mechanics, this degeneracy is lifted through the diagonalization 
of the Hamiltonian, in this case the 'Hamiltonian' is the full Dirac operator. The two "wave 
functions" which diagonalize the "Hamiltonian" are the sum and the difi'erence of the would-be 
zero modes, one of which is a 2-component left-handed spinor, and the other is a 2-component 
right-handed spinor. The resulting wave functions are 4-component Dirac spinors; one can be 
obtained from another by multiplying by the 75 matrix. As the result the two would-be zero 
eigenstates are split symmetrically into two 4-component Dirac states with non-zero eigenvalues 
equal to the overlap integral between the original states: 



0(x-2;l,pl)(a;-2;l)+^7lV^ 

(j){x - Z2, P2){x - Z2)iP^j,e^^, 

P 

7r(2a;2)V2(a;2 + p2)3/2- 



(4.15) 



A = ±|Tii|, 

2piP2, 
^12 



Til = /d4x$l(x-Zl,C/l)tH^)$2(x-^2,C/2)'^-^°°-^'IV(C/|C/2^12;.<J+). (4.16) 
J Zto 



We see that the splitting between the would-be zero modes falls off as the third power of the 
distance between / and /; it also depends on their relative orientation. The fact that two levels 
have eigenvalues ±A is in perfect agreement with the 75 invariance mentioned in the previous 
section. 

When one adds more /'s and Ps each of them brings in a would-be zero mode. After the 
diagonalization they get split symmetrically with respect to the A = axis. Eventually, for an 
// ensemble one gets a continuous band spectrum with a spectral density i/(A) which is even 
in A and finite at A = 0. 

Let the total number of Ps and Ps in the 4-dimensional volume V be A^. The spread k, of 
the band spectrum of the would-be zero modes is given by their average overlap (4.16): 

^ - yjd zudUu \T,2\ -ywjj^, - 6-62107 - (4.17) 
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Figure 5: Schematic eigenvalue distribution of the Dirac operator. The soUd hues are the zero 
mode and free contributions, the dashed hne an estimate of the full spectrum. 



where 



F{kp) = 2t 



Io{t)K,{t) - h{t)Ko{t) - jIi{t)K,{t) 



{kpf 



F(0) = 1, (4.18) 



is the Fourier transform of the instanton zero mode (4.15); the modified Bessel functions are 
involved here. Numerically, if one takes the instanton density N/V — (1 fm)~^ and the average 
instanton size p = | fm (the old Shuryak's values consistent with the variational estimate 
through A, see the previous section) one obtains k = 100 MeV. 

In the random instanton ensemble, one gets the following spectral density of the Dirac 
operator [11]: 



TTK V 4k^ 

From eq. (4.13) one immediately finds the value of the chiral condensate: 



Nr. 



= -(253 MeV)^ 



(4.19) 



(4.20) 



i?2pV 6.62 

which is quite close to the phenomenological value ^. 

I would hke to stress that the chiral consensate is not hnear in the instanton density N/V 
what one would naively expect but rather proportional to its square root (the gluon condensate 
is, naturally, linear). If the instanton density goes to zero the spectral density of the Dirac 
operator tends to a 5-function at zero eigenvalues. This is what one expects from the zero 
modes in the infinitely-dilute limit. 

Eq. (4.19) is known as the Wigner semicircle spectrum. For high eigenvalues A » k the 
spectral density is asymptotically given by that of free massless quarks: 

^Contrary to the gluon condensate, the chiral or quark condensate is somewhat dependent on the scale where 
one estimates it. The above number refers to the scale given by the average instanton size, that is 600 MeV. 
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.(A)«.^A3. (4.21) 

Schematically, the combination of the low- and high-energy spectra are shown in Fig. 5 where 
the intereference with the intermediate modes with \ >2k has been ignored. 

We see thus that the spontaneous chiral symmetry breaking by instantons can be interpreted 
as a delocahzation of the "would-be" zero modes, induced by the background instantons, re- 
sulting from quarks hopping between them [10, 11]. Imagine random impurities (atoms) spread 
over a sample with final density, such that each atom has a localized bound state for an electron. 
Due to the overlap of those localized electron states belonging to individual atoms, the levels 
are split into a band, and the electrons become delocalized. That means conductivity of the 
sample, the so-called Mott-Anderson conductivity. In our case the localized zero quark modes 
of individual instantons randomly spread over the volume get delocalized due to their overlap, 
which means chiral symmetry breaking. 

This analogy between chiral symmetry breaking in QCD and the problem of electrons in 
condensed matter systems with random impurities goes even further [11]. The acquisition of 
a dynamical mass by a quark is fully analogous to the appearance in the Green function of 
an electron in a metal with impurities of a finite relaxation time (but in our case this time 
depends on the momentum) . The appearance of the massless pole in the pseudoscalar channel 
corresponding to the Goldstone pion is analogous to the formation of a diffusion mode in the 
density-density correlation function. For the recent development of these and related ideas, see 
refs. [72, 73, 74] and references therein. 

Recently, the instanton mechanism of the SCSB has been scrutinized by direct lattice meth- 
ods [75, 76, 77]. At present there is one group [77] challenging the instanton mechanism. How- 
ever, the density of alternative 'local structures' found there explodes as the lattice spacing 
decreases, and this must be sorted out first. Studies by other groups [75, 76] support the 
mechanism described above. 

4.3 Quark propagator and dynamical quark mass 

The spectral density of the Dirac operator, averaged over the instanton vacuum, carries very 
limited information, although one can already see that chiral symmetry is spontaneously broken. 
In fact, one can compute analytically much more complicated correlation functions in the 
instanton vacuum, such as the quark propagator and correlators of mesonic currents [11]. What 
is difficult to calculate analytically can be done by numeric simulations of the instanton ensemble 
[19]. 

Each time a quark 'hops' from one random instanton to an anti-instanton (and vice versa) 
it has to change its helicity, because instanton's zero mode is right-handed while the anti- 
instanton's one is left-handed, see the schematic drawing in Fig. 6. Delocahzation imphes 
quarks make an infinite number of such jumps. An infinite number of helicity-fiip transitions 
generates a non-slash term in the quark propagator, i.e. the dynamically-generated mass M{p), 
see Fig. 7. It implies the spontaneous chiral symmetry breaking. 

Mathematically, one has to consider the quark propagator in the gluon background being 
the superposition of an infinite number of Fs and J's , and then average the propagator over 
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Figure 6: Quarks hopping from instantons to anti-instantons and vice versa flip helicity. An 
inflnite number of such jumps generates the dynamical mass M{p). Actually, one has to take 
into account that quarks can 'return' to the same pseudoparticle inflnitely many times [11]. 



the positions, sizes and orientations of instantons according to their partition function (2.50). 
This is a hopeless task, unless one exploits the fact that the packing fraction of instantons is 
small. The actual expansion parameter is a = n'^p^N/iVN^ ~ 1/20 which is not so bad. In 
the leading order in that parameter one can derive a closed equation for the quark propagator 
averaged over the ensemble. Its solution has the form [11, 13] 

The 'wave function renormalization' factor Z{p) differs from unity by a function proportional 
to the above small parameter a, and this difference will be neglected. The dynamical quark 
mass M{p) is, on the contrary, proportional to the square root of the packing fraction: 



M(/) = const. J^^^^F\pp), (4.23) 



VN, 



c 



with the function F{t) given by eq. (4.18); it is related to the Fourier transform of the zero 
mode ^. 

The overall numerical constant is found from the self-consistency or gap equation [11]: 

r d^p M\p) _ N 
"J (27r)4M2(p)+p2 - y- ^4-24) 

For the 'standard' values of the instanton ensemble, N/V = (lfm)~^, p = |fm, one gets at 
zero momentum M(0) = 345 MeV. The dynamical mass (4.23) is plotted in Fig. 7 on top of 
the recent lattice data for this quantity obtained by an extrapolation to the chiral limit [78] . 

Knowing the quark propagator one is able to compute the chiral condensate directly with- 
out referring to the Banks-Casher relation. By deflnition, the chiral condensate is the quark 
propagator taken at one point; in momentum space it is a closed quark loop: 

^It has been known from the perturbative analysis of the 1970's that asymptoticaUy M{p) ~ {as /in) <tpip> 
/p^ whereas eq. (4.23) gives at large virtuality M{p) ~ This is because perturbative gluons are ignored 

in the instanton derivation. At very large p the perturbative regime ^ Ijp^ has to take over. 
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Figure 7: Dynamical quark mass M{p) [11] as function of quark virtuality. The scattered points 
are lattice data obtained by exprapolation to the continuum and chiral limits [78]. Courtesy 
P. Bowman. 



r d*v M(v) I NN 

-(*«M,„. = .{^**^«E., = iN. / = const. (4.25) 

Putting in the 'standard' instanton ensemble parameters one gets the same value of the con- 
densate as before: (ipip) = — (253MeV)^. 

Furthermore, using the small packing fraction as an expansion parameter one can also 
compute [11] more complicated quantities like 2- or 3-point mesonic correlation functions of 
the type 

{JAix)JBiy)), {JAix)JBiy)Jciz)), Ja = ^Ta^Ij (4.26) 

where is a unit matrix in colour but an arbitrary matrix in flavour and spin. Instantons 
influence the correlation functions in two ways: i) the quark and antiquark propagators get 
dressed and obtain the dynamical mass, as in eq. (4.22), ii) quark and antiquark may scatter 
simultaneously on the same pseudoparticle; that leads to certain correlations between quarks 
or, in other words, to effective quark interactions. These interactions are strongly dependent 
on the quark-antiquark quantum numbers: they are strong and attractive in the scalar and 
especially in the pseudoscalar and the axial channels, and rather weak in the vector and tensor 
channels. I shall derive these interactions in the next section, but already now we can discuss 
the pseudoscalar and the axial isovector channels. These are the channels where the pion 
shows up as an intermediate state. Since we have already obtained chiral symmetry breaking 
by studying a single quark propagator in the instanton vacuum, we are doomed to have a 
massless Goldstone pion in the appropriate correlation functions. However, it is instructive to 
follow how does the Goldstone theorem manifest itself in the instanton vacuum. It appears 
that technologically it follows from a kind of detailed balance in the pseudoscalar channel (such 
kind of equations are encountered in perturbative QCD where there is a delicate cancellation 
between real and virtual gluon emission). Since we have a concrete dynamical realization of 
chiral symmetry breaking we can not only check the general Ward identities of the partially 
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conserved axial currents (which work of course) but we are in a position to find quantities 
whose values do not follow from general relations. One of the most important quantities is the 
constant: it can be calculated as the residue of the pion pole. One obtains [11, 79]: 



V (27r)4 [M2(p)+p2]2 
N R 

= const. • —pMn- fa (100 MeV)2 vs. mUeVf (exper.). (4.27) 
V p 

This is an instructive formula. The point is, is anomalously small in the strong inter- 
actions scale which, in the instanton vacuum, is given by the average size of pseudoparticles, 
1/p ~ 600 MeV. The above formula says that F^^ is down by the packing fraction factor 
(p/Ry ~ 1/9. It can be said that Fj^ measures the diluteness of the instanton vacuum. How- 
ever it would be wrong to say that instantons are in a dilute gas phase - the interactions are 
crucial to stabilize the medium and to support the known renormalization properties of the 
theory, therefore they are rather in a liquid phase, however dilute it may turn to be. By calcu- 
lating three-point correlation functions in the instanton vacuum it is possible to determine e.g. 
the charge radius of the pion as the Goldstone excitation [11]: 

~ (340MeV)-^ vs. (310 MeV)"^ (exper.). (4.28) 

In flavour-singlet pseudoscalar channel the instanton-induced interactions are not strong 
attraction as in the non-singlet channel. Therefore, the 1]' meson is not a Goldstone boson: 
the famous Ua{1) problem is solved by instantons, as anticipated at the very beginning of the 
instanton era by 't Hooft [3]. Moreover, in the limit Nf/Nc — >■ instantons reproduce [12, 27] 
the theoretical Witten-Veneziano formula [80, 81] for the singlet rj' mass, as given by 

TT 

where < > /V =< {N^ — N_y > /V is the topological susceptibility. In the instanton 
vacuum it is related to the difference between the number of /'s and /'s . It should be stressed 
that eq. (4.29) is correct in the chiral hmit, and that the topological susceptibihty is that of 
the pure-glue world, without quarks. As emphasized in subsection 2.7, the instanton vacuum 
is described by the grand canonical ensemble of J's and J's , with the fluctuating number of 
pseudoparticles N±. In the CP-conserving vacuum, i.e. for the 'instanton angle' 6* = 0, one 
flnds the equal averages < >—< >~ V from a saddle-point equation. At 7^ the 
saddle-point values for < N± > are complex conjugate to each other [7, 27] . The square of the 
difference between the numbers of Fs and J's is behaving in the normal thermodynamic way, 
<{N^ — iV_)^>~ V, and gives rise to the topological susceptibflity <Q^>. 

When the back influence of quarks on the instanton ensemble is taken into account, which is 
a 0{Nf/Nc) effect, the average < {N^ — N_)'^ > gets dynamically suppressed since at 7^ N_ 
the number of left- and right-handed zero modes are not equal, and the fermion determinant 
goes to zero in the chiral limit [82, 83, 27]. This is in accordance with the general anomalous 
Ward identites in the Ui{A) channel [81]. 
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The Witten-Vcneziano formula (4.29) is an idealization at Nc oo, m — > 0. For non-zero 
quark masses there is a singlet-octet mixing [81, 84] resulting in physical rj, rj' mesons. Actually 
the mixing angle appears to be rather small - about 10°. 

Let me note that all quantities exhibit the natural behaviour in the number of colours 
Nc [12, 68]: 

p = 0(1), 7V/(0) = O(l), y^ = 0(l), etc. (4.30) 

A systematic numerical study of various correlation functions in the instanton vacuum has 
been performed by Shuryak, Verbaarschot and Schafer [19], see also the review [22]. In all cases 
considered the results agree well or very well with experiments and phenomenology. 



5 Instanton-induced interactions 

There are two philosophically different but mathematically equivalent ways of computing ob- 
servables in the instanton vacuum. The first is to compute an observable for a given config- 
uration of /'s and /'s and then average over the ensemble [11]. The second is opposite: first 
average over the ensemble [12]. Since two or more quarks can scatter off the same instanton 
or anti-instanton, averaging over their positions and orientations in colour space induce certain 
correlations between quarks, which one can also call the interactions. The result is an effective 
action for quarks which contains instanton effects in induced multi-quark interactions. The 
would-be zero modes serve here as a bridge, passing information from the instanton vacuum 
to the effective quarks through the induced vertex. The consequent interactions are vertices 
involving 2Nf quarks, commonly cited as 't Hooft interactions since he was the first to specify 
the proper quantum numbers. 

It is convenient to decompose the 4-component Dirac bi-spinors describing quark fields into 
left- and right-handed Weyl spinors which we denote as 

^HRy i'l{R)Jai^ (5-1) 

where / = l...Nf are flavour, a = l...Nc are colour and i = 1,2 are spinor indices. Let us 
introduce the 't Hooft-like 2A^j-fermion vertices generated by /'s and /'s , which we denote 
by Yp^\ respectively. These vertices are obtained by explicit averaging over (anti)instanton 
orientation matrices and over the instanton size distribution d{p) (2.54). Averaging over 
instanton positions in 4d Euclidean space-time produces the overall conservation of momenta 
of quarks entering the vertex Y, hence it is convenient to write down the quark interaction 
vertex in the momentum space. There are formfactor functions F{kp) (4.18) associated with 
the Fourier transform of the fermion zero modes of one instanton, attached to each quark 
line entering the vertex. The 2A^/-fermion vertex induced by an instanton (see Fig. 8) is, in 
momentum space. 
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... 

Figure 8: Instanton-induced 2iV/-quark vertex. The black blobs denote the formfactor functions 
F{kp) attached to each quark leg. 



■{2^)%h + ... + -h- ... - iNf) ■ U^U^/e'^'fe,^,^ l^i^lfaf^^kf)d'''' ih)]] ■ (5.2) 

For the Y~ vertices induced by /'s one has to replace left-handed Weyl spinors V'l by right- 
handed ones, ipR, t/'Ij. Using these vertices one can write down the partition function to which 
QCD is reduced at low momenta, as a functional integral over quark fields [12, 71, 27]: 



Nf 



N+ / \ Ar_ 



^ = /Wexp(/Ai;^,.#/J^^] (^J (5.3) 

where are the number of J's and J's in the whole 4d volume V. The volume factors in 
the denominators arise because of averaging over individual instanton positions, and certain 
mass factors Mi ^ are put in to make eq. (5.3) dimensionless. Actually, the mass parameter Mi 
plays the role of separating high-frequency part of the fermion determinant in the instanton 
background from the low-frequency part considered here. Its concrete value is irrelevant for 
the derivation of the low-energy effective action performed below; in fact it is established from 
smooth matching of high- and low-frequency contributions to the full fermion determinant in 
the instanton vacuum [10, 11]. 

Having fermion interactions in the pre-exponent of the partition function is not convenient: 
one should rather have the interactions in the exponent, together with the kinetic energy term. 
This can be achieved by rewriting eq. (5.3) with the help of additional integration over Lagrange 
multipliers A±: 



X+VMi' J V X-VMi' J 



Nf I 

+ / d'x p ^fi^^^Jf + A+y+^ + X-Y^^ I . (5.4) 

Since N± V oo integration over X± can be performed by the saddle-point method; the 
result is eq. (5.3) we started from. 



36 



As seen from eq. (5.4), X± plays the role of the coupling constant in the many-quark in- 
teractions. It is very important that their strength is not pre-given but is, rather, determined 
self-consistently from the fermion dynamics itself; in particlular, the saddle-point values of A-t 
depend on the phase quarks assume in the instanton vacuum. In the chiral symmetry broken 
phase the values of A±, as determined by a saddle-point equation, appear to be real. 

To get the 2A^/-fermion vertices (5.2) in a closed form one has to explicitly integrate over 
instanton orientations in colour space. For the 2A^/-fermion vertex one has to average over Nf 
pairs of {U,W). In particular, one has: 



dU U^Ul^ 



»2 "ll "22 



etc. 



(5.5) 



In the simplest case of one quark flavour, Nf — 1, the "vertex" (5.2) is just a mass term for 
quarks, 

= jdp^m^pnkp)f[€{k)^p^r{k)\ (5.6) 

with a momentum dependent dynamically-generated mass M{k) given by 

M{k) -^Jdp v{p) [2T:pF{kp)f ^ ^ [27:pF{k-p)f . (5.7) 

In fact, it is exactly the same dynamical mass (4.23) as obtained in refs. [11, 13] by first 
considering the quark propagator in a given configuration of /'s and /'s and then averaging 
over the ensemble. This result comes as follows. 

In order to find the overall scale A of the dynamical mass one has to put (5.6) into eq. (5.4), 
integrate over fermions, and find the minimum of the free energy with respect to A-t. At ^ = 
the QCD vacuum is CP invariant so that = = A^/2 and consequently A+ = A_ = A 
In this case the 75 term in gets cancelled, and the exponent of the partition function (5.4) 
reads: 



-7VlnA + jd'xj^Tr\n\^^ + i^[2T:-pF{kp)] 



= -ArinA + 2Ar,F/^ln|A;2+ (A[2^pF(fcp)Py|. (5.8) 

Differentiating it with respect to A and using eq. (5.7) one gets the gap equation eq. (4.24) 
which is in fact a requirement on the overall scale of the constituent quark mass M{k); its 
momentum dependence is anyhow given by eq. (5.7). Since the momentum integration in 

''Fluctuations of the topological charge, A^_|_ — N_, leading to the topological susceptibility (related to the 
solution of the U{1) problem) has been considered in refs. [12, 27]. 
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eq. (4.24) is well convergent and is actually cut at momenta k ~ the saddle-point value 



of the Lagrange multiplier A is of the order of \J N^N /V / p. The steepness of the saddle-point 
integration in A is guaranteed by the volume, hence the use of the saddle-point method is ab- 
solutely justified. Note that eq. (5.6) reproduces the massive quark propagator (4.22), hence 
the chiral condensate is, as before, given by eq. (4.25). 

For two flavours, Nf = 2 averaging eq. (5.6) over the instanton orientations with the help 
of eq. (5.5) gives a nontrivial 4-fermion interaction. It is, of course, non-local: a formfactor 
function F{kp) is attributed to each fermion entering the vertex; in addition it should be 
averaged over the sizes of instantons. The non-locality is thus of the order of the average 
instanton size in the vacuum. One has [12, 79]: 

= I + k2-h- h) 

■ jdpv{p) i27rp)^F{k,p)F{k2p)F{hp)F{kp) 



SN, 



W bPlfi (^i) ^M^V^f (^i)] bPif. ihMf M] . (5.9) 



For the /-induced vertex Y~ one has to replace left-handed components by right-handed ones. 
In all square brackets summation over colour is understood. Note that the last-line (tensor) 
term is suppressed at large N^, it, however, is crucial at Nc = 2 to support the actual SU{4) 
chiral symmetry in that case [12]. The antisymmetric ^^^^'^^g^g.^ structure demonstrates that the 
interactions have a determinant form in the two flavours. Using the identity 

26/^/^6,,,, = - {rX (ry^l (5.10) 
and adding the /-induced vertex one can rewrite the leading- A^c (first) term of eq. (5.9) as 

(V'V)' + (V'SV')' - (V'V^V)' - (V'V^TsV')' (5.11) 

which resembles closely the interaction of Nambu-Jona-Lasinio [85] model. It should be stressed 
though that in contrast to that at hoc model the full 't Hooft interaction i) violates the Ua{1) 
symmetry in a very definite manner, ii) has a fixed interaction strength related to the density 
of instantons and in) contains an intrinsic ultraviolet cutoff due to the formfactor functions 
F{kp). In addition, at A^c = 2 it correctly preserves the actual Pauli-Giirsey SU{A) chiral 
symmetry [12]. It should be added that the naive addition of a nonzero current quark mass 
to the NJL Lagrangian fails to reproduce several known low-energy Ward identities, as well 
as the phenomenologically-known coefficients in the Gasser-Leutwyler chiral Lagrangian (the 
terms containing and m-p^). The microscopic instanton approach preserving all symmetries 
of QCD is capable to correctly incorporate nonzero quark masses, and it does so in a rather 
nontrivial way [86]. 
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Higher number of flavours, Nf > 3, have been considered in ref. [12, 79] where also the 
bosonization of the instanton-induced interactions has been performed. For three quark flavours 
four-fermion interactions cannot be sufficient, as the Ua{1) symmetry is not broken without a 
6-fermion interaction involving all three flavours. In this case the four-fermion interaction is 
not 'fundamental' but arises from the 2A^/-fermion 't Hooft vertex after one integrates out, say, 
the strange quark. 

I have been often asked: Why mention instantons when using four-fermion interactions? 
Isn't it just the Nambu-Jona-Lasinio model? The answer is that the NJL model is not QCD: if 
it is the correct low-energy effective theory it should be derived from QCD. The model has been 
proposed a decade before the advent of QCD and formulated in terms of nucleons, not quarks, 
and time has come to learn what is QCD suggesting instead of the NJL model. A concrete 
form of NJL-like interactions of quarks, with specific formfactors and strength related to A, is 
derived from QCD via instantons. As indicated above, it preserves correctly the symmetries 
of QCD. Phenomenologically, it is successful. Instantons enable one to understand clearly the 
domain of applicability of the low-energy effective theory, its limitations and ways to generalize 
it, for example if one wishes to take non-zero current quark masses or to include perturbative 
gluons, see below. I can admit that the effective four-fermion interactions with some couplings 
and formfactors may follow from the microscopic QCD via another, non-instanton, formalism. 
Had there been an alternative derivation (from dyons, vortices or whatever) then we could com- 
pare the effective models following from different microscopic considerations and choose the one 
which fits experiments better. We do not even know if potential alternative mechanisms of the 
SCSB could lead to anything resembling the NJL model at low momenta. Today there exists 
only one derivation of the NJL-like model and it is from instantons. Therefore, unless the 
situation changes, I would refer to instantons when four-fermion interactions are used as an 
effective low-energy theory. 

There has been an interesting application of the instanton-induced interactions to the prob- 
lem of a phase transition from the usual chiral-broken phase to the so-called colour supercon- 
ducting phase. The point is, the instanton-induced vertex (5.9), being Fierz-transformed to the 
diquark channel (as contrasted to the quark- antiquark one), gives an attraction for two quarks 
with the L — 0, S = 0, isospin = quantum numbers [87, 88]. At A^'c = 2 the strength of this 
attraction is exactly the same as in the qq channel owing to the SU (4) symmetry mentioned 
above [12]. At A^c ^ 3 the attraction in the qq channel prevails, leading under normal cir- 
cumstances to the qq condensation which breaks chiral symmetry [88] . However, at sufficiently 
high baryon number density the relative strengths of the effective attraction in the qq and qq 
channels reverse: it becomes favourable for the diquarks to condense [89, 90]. Since diquarks 
are in a triplet state for the SU (3) colour, diquark condensation means the spontaneous break- 
ing of colour symmetry, hence the term 'colour superconductivity'. On the contrary, the chiral 
symmetry gets restored at high matter density. It is a first order phase transition. 

The instanton-induced phase transition has been studied in some detail in ref. [91] along 
the lines of the present section. To that end one needs to generalize the formfactor functions of 
the instanton-induced vertex to nonzero chemical potentials. The diquark condensation at high 
matter density mimics the Higgs phenomenon and gives rise to 'Meissner' masses of gluons. To 
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find the gluon masses one needs to construct a conserved Nother current from the non-local 
't Hooft vertex or, in other words, to find out how to generalize 't Hooft vertex to include 
perturbative gluons into it [92]. 



6 Effective Chiral Lagrangian 

It is very important that the dynamical quark mass is parametrically much less than 1/p: the 
dimensionless quantity 

{M-pf - ^ « 1 (6.1) 

is suppressed by the packing fraction of instantons in the vacuum. The whole approach to 
the instanton vacuum implies that instantons are on the average relatively dilute and that this 
packing fraction is numerically small. At low momenta k < 1/p = 600 MeV there are exactly 
two degrees of freedom left: quarks with a dynamical mass M <^ 1/p and the massless Goldstone 
pions. Let us formulate the theory to which QCD is reduced at low momenta k < 1/p: it will 
involve only quarks with dynamical mass, interacting with pions. 

Bosonizing the instanton-induced interaction (5.9) and neglecting all heavy fields one arrives 
to the effective theory to which QCD is reduced at low momenta. The theory is defined by the 
partition function [11, 12] 

r r r ( rd^kd^l i 

■ 'i^Uk) (yiix)^-^ + ui^ix)^-^) r'-ii)] ] , 

Ul{x) = (exp^7^^(a;)A^))^ (6.2) 

Eq. (6.2) shows quarks interacting with chiral fields U{x), with formfactor functions equal to 
the square root of the dynamical quark mass attributed to each vertex where U {x) applies. The 
matrix entering in the parentheses is actually a Nf x Nf matrix in fiavour and a 4 x 4 matrix 
in Dirac indices. It can be identically rewritten as 

U{x)^-^ + U\x)^-^ = exp (i7r^(x)A^5) = U^'{x), (6.3) 

the industrious final abbreviation being due to Pavel Pobylitsa. 

The formfactor functions M{kp) for each quark line attached to the chiral vertex auto- 
matically cut off momenta at /c > 1/p. In the range of quark momenta k <^ 1/p (which we 
shall be mostly interested in) one can neglect this non-locality, and the partition function (6.2) 
is simphfied to a local field theory: 

Z = J Dtt^ J Di/j^ Di/j exp J d^x i/j^ (x) [i^ + iMW' {x)]ilj{x). (6.4) 
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One should remember, however, to cut the quark loop integrals at /c ~ 1/p ~ 600 MeV. Notice 
that there is no kinetic energy term for pions: it appears only after one integrates over the quark 
loop, see below. Summation over colour is assumed in the exponent of eq. (6.4). Eq. (6.4) defines 
a simple and elegant local field theory although it is still a highly non-trivial one. This theory 
has been reviewed in some detail in ref . [79] . 



7 Chromo magnetic pomeron 

The soft part of hadron-hadron collisions at high energies is still full of riddles. Why is the 
additive quark model so successful phenomenologically, stating that hadron-hadron coUision 
can be viewed, to a good approximation, as an incoherent sum of collisions of their constituent 
quarks [93, 94]? Why is the total pp cross section around 40 mb? This number, as any other in 
strong interactions, needs to be obtained from the transmutation of dimensions and expressed 
through A (1.1). Why is the differential elastic cross section approximately proportional to 
the fourth power of the electromagnetic formf actor? Why are spin effects represented, e.g., by 
the analyzing power and the polarization of the inclusively produced hyperons, so strong and 
independent of energy [95]? 

The traditional view on a pomeron is that it is basically a two-gluon exchange a la Low- 
Nussinov [96] or its more sophisticated Balitsky-Fadin-Kuraev-Lipatov sibling [97], see the 
recent book [98] for a review. Both versions are infrared-unsafe, in the sense that the transverse 
momenta carried by gluons are cut from below not "by themselves" but rather by the inverse 
sizes of the colliding hadrons. [In the BFKL case one eventualy runs into the infrared region even 
if the hadron sizes are vanishing.] In the Nachtmann-Dosch model [99] where non-perturbative 
scattering of strings between quarks inside hadrons is considered, the characteristic transverse 
scale is also set by hadron sizes; both models are opposite in spirit to the additive quark model, 
so that e.g. the famous Levin-Prankfurt ratio [93] a^N '■ cTt^n '■ crnn = 9:6:4 appears to be 
accidental. None of the pomeron models suggested so far are able to explain coherently the 
riddles of large spin effects in high energy production [100]. 

In this section we would like to suggest another point of view on the soft pomeron From 
the instantons' point of view, the leading soft gluon exchange flips quark's helicity and not 
conserves it, as in perturbation theory. An interference between the two processes may lead to 
the interesting zoo observed in spin physics in high-energy hadron reactions [100]. 

We know already from the previous sections that instantons induce very strong helicity-flip 
interactions of quarks, known as 't Hooft interactions. Being point-like (more precisely, of 
the range of p), they are irrelevant at high energies. However, since 't Hooft interactions are 
induced by instantons that are lumps of gluon field, there exist more complicated helicity-flip 
processes, with the classical gluon field of instantons radiated from quark vertices, see Fig. 9. 
That kind of amplitudes have been studied in relation to baryon number violating processes in 
electroweak theory [101, 31]. 

The simplest helicity-fiip amplitude induced by instantons is presented in Fig. 9a. Its first 

®This section has been prepared in co-authorship with Maxim Polyakov. We thank V. Petrov, M. Ryskin 
and M. Strikman for discussing these matters. 
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b) c) 



Figure 9: Helicity-flip quark-gluon vertices. 

term in the expansion in the gluon field is the quark anomalous chromomagnetic moment first 
introduced by Kochelev [102]: 



^magn = ^ (V; a,^ ^P) G%, a,^ = - [7,7.] , (7. 1) 



where is the gluon field strength, ip is the quark field, M = M(0) is the dynamical quark 
mass at zero virtuality and fi is the anomalous chromomagnetic moment induced by instantons. 



^, = -l^-liM-pr^^^-0.7AA. (7.2) 
2 as — 1 

Although we agree qualitatively with ref. [102], our formula for /i is different and the 
value is several times bigger because of a too crude estimate used there. In our estimates we 
use the formalism and the instanton parameters of refs. [5, 11, 27]: M ^ 345MeV, ~ 
600 MeV, 271 /as ~ 12, n ^ (200 MeV)"^, where n = N/V is the instanton density. We note that 
fi is stable at large number of colours, which is also the property of the perturbation theory. 
The fact that /i is large is not accidental but follows from the relation 

o /1 27r 1 

being the condition for the instanton ensemble to stabilize itself [27]. 

The chromomagnetic vertex (7.1) can be obtained from more general instanton- induced 
vertices of Fig. 9b, c by integrating out quark fields. The general vertex in case of A'^^ fiavours 
is 

+ (75^-75) (7.4) 

where denotes the (fg) minor of the fiavour matrix = ipf^^^ipg. The summation 
over fiavour indices f,g and colour is understood. After bosonization at low momenta [12, 79] 
the general 2A^j-fermion vertex with an additional gluon (7.4) may be presented in a multi- 
Goldstone form [103]: 



9}^ J. 
2M 

In this form the vertex is explicitly chirally invariant. 
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Figure 10: Leading contributions to the inelastic cross section at intermediate energies. 

The large anomalous chromomagnetic moment gives rise to a sizable (possibly dominant) 
contribution to the hadron-hadron scattering at high energies. We would like to stress that 
the quark-gluon vertex (7.1) is gauge-invariant - in sharp contrast to the usual perturbative 
vertex tp'-^^A^ip. To get a gauge-invariant result, that vertex has to be inserted in all possible 
quark lines belonging to a 'colourless' hadron. In other words, the 7^ vertex interacts with the 
quarks' colour charge, and one must take into account that hadrons have zero colour charge. 
That is why the Low-Nussinov-BFKL pomeron is sensitive to the size of hadrons. For the 
chromomagnetic vertex, the colour interference is not inevitable. Moreover, it is suppressed 
as ~ 1/ {k±R) where R is the hadron size and k±_ is the typical momentum transfer through 
instantons, which is about 1 GeV, as we shall see. Therefore, one can neglect, to the zero 
approximation in this small parameter, the interference of diagrams where chromomagnetic 
vertices are attached to different quarks in the amplitude and in the conjugate amplitude. 
This immediately leads to the conclusion that hadron cross sections are incoherent sums of the 
quark-quark cross sections, and that is the basis of the additive quark model. 

Let us estimate the imaginary part of the forward scattering amplitude corresponding to the 
simplest diagram of Fig. 10a. The gluon exchange in the t-channel is in fact a complicated ob- 
ject as it corresponds to the interaction of two coherent classical gluon fields of (anti) instantons 
in a non-linear Yang-Mills theory. The stringy line in Fig. 10 is actually the connected part 
of the correlation function < G^^(x)G'^^(y) > where both field strengths are superpositions of 
(anti) instantons belonging to the two vertices. 

Assuming the sum Ansatz for the two instanton fields the correlation function can be directly 
calculated at all k using formulas from ref. [7]. A cleaner way to calculate it would be to use 
the unitarity and analiticity method of refs. [31, 32] but that has not been done. For an 
estimate, we use the simplest model, replacing each field strength just by the derivative of 
the instanton potential. The model can be conveniently written in the form of the anomalous 
chromomagnetic vertex supplemented by the instanton formfactor: 



Magnetic vertex 



M 



(7.6) 




(7.7) 
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where K{z) are the modified Bessel functions, A(0) — 1. The fact that there must be some 
formfactor is model-independent as one cannot transfer momenta larger than 1/p to a lump of 
coherent field of size p without destroying it. 

Using the optical theorem, cxtot = ImA(s,t = 0)/s, and neglecting the constituent quark 
masses at large energies s we estimate the total qq cross section arising from Fig. 10a as ^ 

where the last factor comes from contracting the spin traces over unpolarized quarks, 

iTr [^o-^^/l:^.(^ + ^)aaf3kf3] ■ |Tr [:^2(^^,^k^{:^2 - ^)(^axkx], (7.9) 

and using s ^ M^, k^ and k^ ~ — A;^. We see that the cross section is constant in energy which 
is of course due to the spin 1 exchange. Despite the additional factors of ka in magnetic vertices 
as compared to the ^y^ ones, integration over k± converges at ~ 1/p owing to the instanton 
formfactor which we have modelled by the instanton field itself. 

Eq. (7.8) should be multiplied by the colour factor 4(-j^'Li) arising from contracting twice the 
colour matrices t"' in a colourless meson or baryon. It wiil then become the 'true' quark-quark 
cross section which one can multiply by 4, 2Nc or A'^ to get the meson-meson, meson-baryon 
or baryon-baryon cross sections, respectively. 

Integrating over k± we obtain 

The last (numerical) factor is the dimensionless integral over the specific instanton formfactor 
(7.7); it is of the order of unity but its concrete value may change in a more precise theory. 
Also, one has to add the multi-quark processes of Fig. 9b,c leading to Fig. 10b and the kind. 
Although the cross section of Fig. 10b is A'c = 3 times suppressed as compared to that of 
Fig. 10a, a 30%-correction squared moves the estimate (7.10) close to the qq cross section of 
4 mb needed to explain the 36 mb of the inelastic pp and pp cross sections, before they start to 
rise. 

In a wide range of ^/s ~ 5 — 50 GcV where the total cross sections are approximately 
constant and the rapidity plateau is not fully developed, the inelastic collisions are presumably 
dominated by the decays of hadrons whose constituents experience the qq scattering shown 
in Fig. 10. At higher energies, the rapidity plateau develops, possibly due to the instanton- 
ladder mechanism of ref. [104], and that is accompanied by a weak (logarithmic or small-power) 
increase of the cross section. As shown in ref. [32], the point-like instanton- induced cross section 
for multi-gluon production rises at low and decreases at high invariant masses; the interpolation 
of the two curves inside the intermediate mass range gives the position and the height of the 
maximum suspiciously close to what corresponds to the production of the spherically symmetric 
'sphaleron' configuration of gluons with the invariant mass mgph such that [32] 

^More precisely, the cut in Fig. 10a corresponds to the inelastic cross section. 
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2^msphP=g, (7.11) 

from where one can estimate mgph ~ 2.7 GeV and predict its decay spectrum. A possible 
relevance of sphaleron production in high-energy heavy ion collisions has been stressed by 
Shuryak [105]. However, it is not evident that the same sphaleron production survives in the 
ladder kinematics in the pp collisions as assumed in ref. [104], because of the possible gluon 
interference which has not been fully investigated. Whatever the mechanism of the plateau 
formation at very high energies, it gives a relatively weak rise of the total cross section - on 
top of what? This section attempts to answer this basic question. 

To get the elastic hadron cross section or, more generally, the diffractive dissociation, one 
needs to square the diagrams of Fig. 10 multiplying it by the probability that one or both collid- 
ing hadrons collect back the constituent quarks that have been struck. Since the size of hadrons 
is much larger than the range p of the constituent quarks' interaction, one can roughly consider 
it as point-like. Therefore, the probability that a hadron remains one hadron after a qq collision 
is basically given by the same formfactor as in a true point-like electromagnetic collision. Wc 
get, thus, a natural explanation why the differential elastic cross section is proportional to the 
fourth power of the electromagnetic formfactor [94]. Deviation from this remarkable rule is 
expected at momentum transfers greater than ~ 0.6 GeV, and indeed the elastic cross section 
is decreasing faster than at large t. 

We now briefly discuss the ensueing spin effects. Despite the unusual spin structure of the 
gluon exchange the dominant structure of the pomeron as a whole has the usual ® /^i form. 
However, there is always an interference between helicity-flip and non-flip amplitudes. A large 
energy-independent non-flip amplitude is, again, induced by instantons, this time by the part 
of the fermion propagator in the instanton background that does not contain chiral zero modes 
[106]. Estimating the corresponding non-flip vertex on the quark 'mass shell', ~ M^, wc get 
~ g{j^t"') ■ in'^np'^ {27T / as) / which is of the same order as the perturbative vertex owing to the 
estimate (7.3) but effectively {Mpf pa 1/3 of the magnetic vertex (7.1). It is also accompanied 
by some formfactor B{kp). 

The enhancement of spin-ffip amplitudes as compared to the non-ffip ones is of a general 
nature: The former are related to the spontaneous symmetry breaking and arc therefore non- 
analytic in the (small) instanton density, namely they are proportional to the square root of 
the density N/V [11], whereas the latter, being 'normal', are naturally linear in the density. 
It suggests that at small momentum transfer the spin-ffip chromomagnetic exchange may be 
quite large, if not dominant. The above estimate of the total cross section demonstrates that 
it can well explain the bulk of the total cross section. 

At the same time, the non-ffip is not negligible, and the interference leads to sizable and 
energy- independent polarization effects. Speaking generally, we see that as much as a {MpY fa 
30% polarization of quarks can be expected in experiments with unpolarized hadrons. 

As correctly stressed by Soffer [100], the polarization is a purely quantum-mechanical effect: 
it is nonzero only if the relative phase of the spin-flip and non-flip amplitudes is nonzero. 
Therefore, one has to go beyond the Born approximation of Fig. 10. The first diagram leading 
to quark polarization is with three gluon exchanges. If one takes one of the three vertices to 
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be 'electric' (i.e. spin non-flip) and two 'magnetic' (i.e. spin-flip) such a diagram will represent 
an interference between a real (Born) spin-flip amplitude and almost purely imaginary (box) 
non-flip amplitude, or vice versa, so that the relative phase is maximal possible. 

We estimate the quark polarization arising from three-gluon exchange diagrams to be 

P^7r\p^\pCip^p), (7.12) 

where p± is the transverse momentum of the quark whose polarization is studied and C is a 
combination of instanton formfactors A, B and their integrals. At present we have only a model- 
dependent C{p±p) which is not very useful. However, the message is clear: The polarization 
rises linearly in p±, is independent of energy and at p_L ~ 0.6 GeV reaches its maximal value of 
more than 10%. Qualitatively, it seems to describe correctly the main striking features of the 
polarization experiments [95]. 

If a quark of new flavour (not present in the incoming hadron) is produced, like in the 
pp ^ A polarization experiment, it goes via the many- flavour version of the chromomagnetic 
coupling shown in Fig. 9b, c. In the bosonized language of eq. (7.5), it means that one must 
insert twice the 75 matrix in the quark line, which inverses the sign of the polarization with 
respect to the case when flavour is not changed! This observation may explain why A po- 
larizations in K~p and pp production have opposite sign, and possibly other mysteries of the 
polarization / analysing power zoo [95, 100]. 

Finally, it is natural to assume that the same process that dominates soft coUisions is 
responsible for the gluon distribution in hadrons, at least at low virtuality Q 0.6 GeV where 
the hard gluon bremsstrahlung has not yet set in. Substituting the standard DGLAP calculation 
by that with the anomalous magnetic vertex (7.6) we obtain the gluon distribution in the 
constituent quark as function of Bjorken's x: 

' 87r2 2N^ \M) xJ ^ ^ y,/iZr^J 

= ^—^ p. 2 ^1.5^—^, (7.13) 

X X 

the polarized gluon distribution Ag{x) being nonsingular at small x and additionally suppressed 
as {Mpy. Taken hterally, eq. (7.13) says that 75% of the nucleon momentum is carried by glue, 
which is a factor of 2 too much, but that factor arises from the model instanton formfactor used 
here and may be an artifact. The x-dependence of the gluon distrubution is more interesting. 
We get the 1/x behaviour at small x, which is consistent with phenomenological fits at low 
virtuality. 

To conclude this section: Knowing the microscopic mechanism of the SCSB one can reveal 
the non-perturbative glue inside the constituent quark. The largest coupling of the 'soft' glue 
to the massive constituent quark appears to be via its anomalous chromomagnetic moment. 
Therefore, the magnetic gluon exchange is expected to be large if not dominant in high energy 
quark scattering. It seems that this picture is supported by many facts we know about high 
energy collisions. 
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Figure 11: Density-density correlation function in the nucleon [21]. Filled circles are measure- 
ments in the full gluon vacuum (corresponding to Fig. 3a,b) while open circles are measured in 
the vacuum with instantons only (Fig. 3c,d). Despite that linear confining potential is absent 
in the instanton vacuum the nucleon structure seems to be very well reproduced. 

I would like to remark that very recently an interesting suggestion on the role of instantons in 
the transition from hard to soft physics at small x has been made in refs. [107, 108]. Instantons 
can induce specific events in deep inelastic scattering [109] which have been recently studied 
experimentally [110]. 

8 Baryons 

There is a remarkable evidence of the importance of instantons for the baryon structure. In 
ref. [21] the so-called density-density correlation function inside the nucleon has been measured 
on the lattice both in the full vacuum and in the instanton vacuum resulting from the full one 
by means of cooling. The correlation in question is between the densities of u and d quarks 
separated by a distance x inside the nucleon which is created at some time and annihilated at 
a later time. The two correlators ('full' and 'instanton') are depicted in Fig. 11: one observes 
a remarkable agreement between the two, up to a; = 1.7 fm. 

It must be stressed that neither the one-gluon exchange nor the linear confining potential 
present in the full gluon vacuum survive the smoothing of the gluon field shown in the lower 
part of Fig. 3. Nevertheless, quark correlations in the nucleon remain basically unaltered! It 
means that neither the one-gluon exchange nor the linear confining potential are important for 
the quark binding inside the nucleon. As a matter of fact, the same remark can be addressed to 
the lightest mesons vr and p since the density-density correlators for these hadrons also remain 
basically unchanged as one goes from the full glue to the reduced instanton vacuum [21]. 
Therefore, one must be able to explain at least the hghtest tt, p, N on the basis of instantons 
only. 

The dynamics remaining in the instanton vacuum is the SCSB, the appearance of the 
dynamical quark mass M{p), and quark interactions induced by the possibility that they scatter 
off the same instanton. Actualy these interactions named after 't Hooft, are quite strong, see 
section 5. They are in fact so strong that for quark and antiquark in the pion channel they 
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Figure 12: 't Hooft interactions in the nucleon (left) essentially come to quarks interacting via 
pion fields (right). 



eat up the 700 MeV of twice the constituent quark mass to nil, as required by the Goldstone 
theorem. In the vector meson channel 't Hooft interactions are suppressed, and that is why 
the p mass is roughly twice the constituent quark mass. In the nucleon they are fully at work 
but in a rather peculiar way: instanton-induced interactions can be iterated as many times as 
one wishes in the exchanges between quarks, see Fig. 12, left. It can be easily verified that the 
diagram in Fig. 12, left, can be drawn as three continuous quark lines going from the l.h.s of 
the diagram to its r.h.s., without adding closed loops. Therefore, that kind of interaction arises 
already in the so-called quenched approximation. At the same time, it yields plenty of Z-graphs 
absent in "valence QCD" but which are necessary to reproduce hadron properties [111]. 

Summing up all interactions of the kind shown in Fig. 12, left, seems to be a hopeless task. 
Nevertheless, the nucleon binding problem can be solved exactly when two simplifications are 
used. The first exploits the fact that in the instanton vacuum there are two lightest degrees 
of freedom: pions (since they are the Goldstone bosons) and quarks with the dynamical mass 
M. All the rest collective excitations of the instanton vacuum are much heavier, and one may 
wish to neglect them. Pions arise from summing up the qq bubbles schematically shown in Fig. 
12, left. The resulting effective low-energy theory takes the form of the non-linear a-model 
introduced in section 6: 



eff 



(8.1) 



The absence of the explicit kinetic energy term for pions (which would lead to the double 
counting) distinguishes it from the Manohar-Georgi model [112]. Expanding the exponent to 
the first power in we find that the dimensionless pion-constituent quark coupling. 
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M 



(8.2) 



is quite strong. The domain of applicability of the low-energy effective theory (8.1) is restricted 
by momenta k < 1/p = 600 MeV, which is the inverse size of constituent quarks. At higher 
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momenta one starts to feel the internal structure of constituent quarks, and the two lightest 
degrees of freedom of eq. (8.1) become insufficient. However, the expected typical momenta 
of quarks in the nucleon are of the order of M ^ 345 MeV, which is inside the domain of 
applicability of the low-momentum effective theory. 

The chiral interactions of constituent quarks in the nucleon, following from the effective 
theory (8.1), are schematically shown in Fig. 12, right, where quarks are denoted by lines 
with arrows. Notice that, since there is no explicit kinetic energy for pions in eq. (8.1), the 
pion propagates only through quark loops. Quark loops induce also many-quark interactions 
indicated in Fig. 12 as well. We see that the emerging picture is rather far from a simple 
one-pion exchange between the constituent quarks: the non-linear effects in the pion field are 
essential. 

The second simplification is achieved in the limit of large Nc- For A^^ colours the number of 
constituent quarks in a baryon is and all quark loop contributions are also proportional to 
Nc- Therefore, at large Nc one can speak about a classical self- consistent pion field inside the 
nucleon: quantum fluctuations about the classical fleld are suppressed as The problem 

of summing up all diagrams of the type shown in Fig. 12 is thus reduced to flnding a classical 
pion fleld pulling Nc massive quarks together to form a bound state. 

9 Chiral Quark-Soliton Model 

Let us imagine a classical time-independent pion field which is strong and spatially wide enough 
to form a bound-state level in the Dirac equation following from eq. (8.1). The background 
chiral field is colour-neutral, so one can put Nc quarks on the same level in an antisymmetric 
state in colour, i.e. in a colour-singlet state. Thus we obtain a baryon state, as compared to 
the vacuum. 

One has to pay for the creation of this trial pion field, however. Since there are no terms 
depending directly on the pion field in the low-momentum theory (8.1) the energy of the 
pion field is actually encoded in the shift of the lower negative-energy Dirac sea of quarks, as 
compared to the free case with zero pion field. The baryon mass is the sum of the bound-state 
energy and of the aggregate energy of the lower Dirac sea, see Fig. 13. It is a functional of the 
trial pion field; one has to minimize it with respect to that field to find the self-consistent pion 
field that binds quarks inside a baryon. It is a clean-cut problem, and can be solved numerically 
or, approximately, analytically. The description of baryons based on this construction has been 
named the Chiral Quark-Soliton Model (CQSM) [113, 114, 115]. 

The model reminds the large-Z Thomas-Fermi atom where Nc plays the role of Z. Fortu- 
nately, corrections to the model go as l/Nc or even as 1/iV^ and have been computed for many 
observables. In the Thomas-Fermi model of atoms corrections to the self-consistent (electric) 
field are of the order of 1/ y/Z and for that reason are large unless atoms are very heavy. 

In the end of the 80's and the beginning of the 90's dozens of baryon characteristics have been 
computed in the CQSM, including masses, magnetic moments, axial constants, formfactors, 
splittings inside the mutliplets and between multiplcts, polarizability, fraction of nucleon spin 
carried by quarks, etc. - see [116, 79, 117, 118] for reviews and references therein. 

Starting from 1996 a new class of problems have been addressed, namely parton distributions 
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Figure 13: The nucleon mass is a sum of the energy of three 'valence' quarks -Bievei and of 
the infinite number of 'sea' quarks whose energy is the aggregate energy of the lower Dirac 
continuum, i^fieid- Solid lines are occupied states in the pion background field, each Nc times 
degenerate. 

in the nucleon at low virtuality [119]. Parton distributions are a snapshot of the nucleon in 
the infinite momentum frame. One needs an inherently relativistic model in order to describe 
them consistently. For example, a bag model or any other nonrelativistic model with three 
quarks in a bound state, being naively boosted to the infinite-momentum frame gives a negative 
distribution of antiquarks, which is nonsense. On the contrary, being a relativistic field-theoretic 
model CQSM predicts parton distributions that satisfy all general requirements known in full 
QCD, like positivity and sum rules constraints. 

Numerous parton distributions have been computed in the CQSM, mainly by the Bochum 
group [119, 120, 121], see also [122, 123, 124] and references therein. There have been a number 
of mysteries from naive quark models' point of view: the large number of antiquarks already 
at a low virtuality, the 'spin crisis' [125], the large flavour asymmetry of an^zquarks, etc. The 
CQSM explains all those 'mysteries' in a natural way as it incorporates, together with va- 
lence quarks bound by the isospin-1 pion field, the negative-energy Dirac sea. Furthermore, 
the CQSM predicts nontrivial phenomena that have not been observed so far: large fiavour 
asymmetry of the polarized antiquarks [121], transversity distributions [126], peculiar shapes 
of the so-called skewed parton distributions [127] and other phenomena in hard exclusive reac- 
tions [128]. Baryon dynamics is rich and far from naive "three quarks" expectations. 

Finally, I would like to stress that physical quantities discussed in this review (and many 
more remaining beyond) are all expressed, via the 'transmutation of dimensions', through the 
only dimensionful quantity there is in QCD, that is A, as it should be in the ultimate theory. 
Instantons serve as a bridge from the first principles of QCD - via the microscopic mechanism 
of spontaneous chiral symmetry breaking - to the observables. Once the basic properties of the 
instanton vacuum are established one does not need to play with parameter fitting, going from 
one observable to another. Instantons provide a coherent picture of very different phenomena 
in strong interactions, and I find it quite remarkable. 
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10 Summary 

1. The would-be linear confining potential of the pure glue world is necessarily screened by 
pion production at very moderate separations between quarks. Therefore, light hadrons need 
not be sensitive to confinement forces but rather to the dynamics of the spontaneous chiral 
symmetry breaking (SCSB). 

2. Very likely, the SCSB is driven by instantons - large non-perturbative fluctuations of the 
gluon field having the meaning of tunneling. The SCSB is due to 'hopping' of quarks from one 
randomly situated instanton to another, each time flipping the helicity. The instanton theory 
of the SCSB is in agreement with the low-energy phenomenology (c/. the chiral condensate 
< >, the dynamical quark mass M{p), F^, rrir,'...) and seems to be confirmed by direct 
lattice methods. Furthermore, lattice simulations indicate that instantons alone are responsible 
for the properties of lightest hadrons tt, p,N, .... 

3. Instantons induce not only very strong non-perturbative quark interactions but also 
new and interesting vertices with an additional gluon emission. In particular, they induce a 
large anomalous chromomagnetic moment which can play an important role in soft high-energy 
hadron scattering, e.g. in spin phenomena. 

4. Summing up instanton- induced quark interactions in baryons leads to the Chiral Quark- 
Sohton Model where baryons appear to be bound states of constituent quarks pulled together 
by the chiral field. The model enables one to compute numerous parton distributions, as well 
as 'static' characteristics of baryons - with no fitting parameters whatsoever. 

5. For highly excited baryons (m = 1.5 — 3 GeV) the relative importance of confining forces 
vs. those of the SCSB may be reversed. One can view a large-spin J resonance as due to a 
short-time stretch of an unstable string or, alternatively, as a rotating elongated pion cloud 
[18]. What picture is more adcqiiatc is a question to experiment. In the first case the dominant 
decay is on the average of the type Barj Bar^j/2 + Mes^j/2; in the second case it is mainly 
a cascade Barj — > Barj_i -|- tt — > Barj_2 -|- tttt — > ... Studying resonances can elucidate the 
relation between chiral and confining forces. 
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